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"^sf. Abstract 

Consider the capillary water waves equations, set in the whole space with infinite depth, and 
CN ■ consider small data (i.e. sufficiently close to zero velocity, and constant height of the water). We 

prove global existence and scattering. The proof combines in a novel way the energy method with 
a cascade of energy estimates, the space-time resonance method and commuting vector fields. 
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In this manuscript we consider the global existence and asymptotic behavior of surface waves 
for an irrotational, incompressible, and inviscid fluid in the presence of surface tension. The fluid 
velocity is given by Euler's equation in a domain U: 

U = \J&t, 9 = &t = {(x,z) = (x!,x 2 ,z) G M 3 , z < h{x,t)}, 
t 

and the free boundary of the fluid at time t 

SB = SB t = {( x ,h(x,t)),x G K 2 } = d@ 

moves by the normal velocity of the fluid. The surface tension is assumed to be proportional (by 
the coefficient c) to the mean curvature k of 33 and we neglect the presence of gravity. In this 
\& " setting the Euler equation for the fluid velocity v, and the boundary conditions are given by 

t-h ; 

(lla) (D t v= f d t v + V v v = -Vp {x,z)£@, 

(1.1b) 



X 



I d t h + V v (h - z) = (x,z)€@, 
1 p = ck, (x, z) G 3%. 

where V = (d±, 82, d z ). Since the flow is assumed to be irrotational, the Euler equation can be 
reduced to an equation on the boundary and thus the system of equations (E-BC) reduces to a 
system defined on SB. This is achieved by introducing the potential ipn where v = Vipu. Denoting 
the trace of the potential on the free boundary by tp(x, t) = ip%(x, h(x, t),t), the system of equations 
for ip and h are (see for instance |30| ) 

( d t h = G(h)^ 

(WW) <^ dtij, = ck- \\d^ + 2{1+ \ dh?) (G(h)i> + 8h • d^f 

{ M)(t = 0) = (ho,ih). 

where 



G(h)i/> = v 7 ! + \dh\ 2 M{h)il) ~ Aip + (quadratic), 
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M being the Dirichlet-Neumann operator associated with Ql; 8 = (81,82); A = \8\; and where the 
mean curvature can be expressed as 

1 ^ ( 8h \ 1 A , . , . . 
k = —8 • I ; ~ -Ah + (cubic). 

In the sequel we take c = 2 for simplicity. The dispersive nature of (IWWft is revealed by writing 
the linearization of this system around (h,ift) = (0,0): 

n 9 ^ / 8 t h = Aip, 

1 ' \ 8^ = Ah, 

where 9^ are at least quadratic in (h,tp). By setting 

u d = A l/2 h + iij) and u d = A l/2 h + ii[> , 
the above system can be written as a single equation 

8 t u = -iA 3 / 2 u + 9t 



(1.3) 

where SH = A39K X + m z 



u(t = 0) = Uq. 



1.1. Known results. A great number of works has been devoted to local well-posedness of the 
water wave problem. These works consider the problem with gravity, capillarity, or both; it is 
furthermore possible to allow domains of finite depth. The small data problem was first addressed 
by Nalimov [26] (see also H. Yoshihara |35|). but the first breakthrough in solving the local well- 
posedness problem with general data is due to Wu [311 132j . There are many other works on local 
well-posedness: we mention in particular Craig [13], Christodoulou and Lindblad [8], Lannes [25] . 
Coutand and Shkoller [12] , Ambrose and Masmoudi [SHU [3], Shatah and Zeng [29| . Alzazard, Burq 
and Zuily [I]. 

In connection with the local regularity problem, Christianson, Hur and Staffilani [7] and Alazard, 
Burq and Zuily [2] have been able to prove recently a nonlinear smoothing effect. 

Fewer works address the global evolution problem; all results are then restricted to small data. 
The first progress in this direction is due to Wu [33] , who proved almost global existence of gravity 
waves in dimension 2. The authors of the present paper [20j could then prove global existence of 
gravity water waves in dimension 3. A very similar result was obtained later by Wu [34] . using a 
different proof. We shall come back later to the methods of proof employed. 

Another line of research has been concerned with traveling waves. We focus on fully localized 
traveling waves in dimension 3 (that is, two-dimensional interface). Such waves have been con- 
structed recently by Groves and Sun [21], and Buffoni, Groves, Sun and Wahlen [6]. The setting 
is that of gravity-capillary waves of finite depth: denote respectively c, g and h for the capillarity 
coefficient, gravity, and the depth of the fluid, and define the Bond number f3 = The aforemen- 
tioned works show that if (3 > |, traveling waves of arbitrarily small size in 1? exist. This should 
be contrasted with our main theorem here: we prove that for c > 0, g = 0, h = 00, small data lead 
to scattering for large time. In our main theorem however, smallness is expressed in weighted L 2 
spaces, and it is not clear to which weighted L 2 spaces the traveling waves of Buffoni, Groves, Sun 
and Wahlen belong. For a more general account of traveling waves, we refer to the recent book of 
Constantin 1111. 



1.2. Main theorem. To state our main result we need to introduce the following notation: let 
8 d = f {8i,d 2 ), n d = f x 1 d 2 -x 2 8 l = Jdi , and S d = f |tft + x% , 
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and let T denote any of the fields T = S,f2, or <9 3 , where <9 3 = d^dtf with i\ + 12 = 3. We adopt 
the multiindex notation: if 7 = (71,72,73) £ N 3 , 

pi _ £71 Q12Q313 _ 



We shall also need the spatial part of 5 



E d = x i d i . 



I ll d ^ 
\ u \\w k <p 



Hi 1 ; : the classical Sobolev space W k,p , whose norm reads 

1 + A k ^j u 

For k a positive real number and £ an integer, define the weighted Sobolev space W^' p by its norm 

Hlvr** = X H r7u lliy fe .p • 

l7l<< 

Recall that u = A x l 2 h + iip, where (h, ip) are given by ((WW]). 
Theorem 1.1. Assume that the initial data uq satisfies 
( L4 ) ll^ 1/2u o|| w 9/2,2, TO2 , + ||yr*«o||wO.?/ M 2^ < £ > 



where K > 10, and e, a* > are sufficiently small. Then there exists a global solution u of (|1.3p 
such that \\u\\x % £ (where the \\ \\x norm is defined in section WJfy . Furthermore, this solution 
scatters, i.e. there exists a solution ui of the linearized problem 

d t U£ = -iA 3/2 u e 

such that 



A 1 ' 2 (u(t) - Ul {t)) 



— > as t — > 00. 

2 



Remark 1.1. Consider ho £ 5(M 2 ) and Vq £ <S(T>q) such that cwclv = divv = 0. Then the data 
fi(ho,vo) satisfy the hypotheses of the theorem if fi is small enough. However note that if we took 
a* = in the theorem above, namely that H^oH^/O^^) < e then n(ho,Vo) would not satisfy the 
hypotheses of the theorem unless a moment condition is satisfied by Vq. Thus by requiring a* > 
we allow a larger class of data at the expense of complicating the proof, since such a condition on 
the initial data makes dh decay at a polynomial rate slower than 1/t (even at the linear level!). 

The proof of this result is based on combining in a new way 1) the vector field method (see 
Klainerman [23]) which is based on invariances of the equation, 2) A cascade of energy estimates, 
3) dispersive bounds and 4) the space time resonance method [181 120) ■ which is based on resonant 
interactions of waves. Below we give a brief and simplified description of the ideas of the proof. 

1.3. Discussion of the difficulties and the method. There are several difficulties that distin- 
guish this problem from the gravity water waves problem, namely: 

1) Low frequencies. The group velocity for capillary waves is given by flCI 1 ^ 2 while for gravity 
water waves it is given by ^|£l -1 ^ 2 - This makes high frequency capillary waves disperse faster 
than gravity waves, but it also causes low frequency waves to decay very slowly. For the long 
time behavior the low rate of dispersion of low frequencies is more problematic. In particular, in 
our setting, the quantity ||V/i||x,oo decays at a polynomial rate slower than 1/t. This slow decay 
combined with the commutation properties of the vector fields causes a cascading growth on the 
weighted energy estimates. 

2) Quadratic resonances. A second difficulty is linked with the fact that £ — > |£| 3//2 is convex rather 
than concave as in the gravity water waves problem. Hence there are non trivial (time) quadratic 
resonances. This was not the case for the gravity water waves. 
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3) Weighted energy estimates. In this work, we chose to estimate u in weighted energy spaces given 
by the invariances of the equation, rather than estimating, say xf in L 2 . This gives better control, 
since the invariant vector fields have better commuting properties than commuting x in the equation 
of the profile. The presence of quadratic time resonances in this problem also makes bounding xf 
in I? a delicate question. 

Through the use of vector fields, we actually integrate the vector fields method and the space 
time resonance method to prove our result. The steps in our proof are reminiscent of the proof of 
global existence in three dimensions for nonlinear wave equations satisfying the null condition by 
Klainerman [23], see also Sogge [27]. In these works, the null condition is utilized via a pointwise 
bound for the quadratic term of the type (we denote Q for the quadratic term, and r for the vector 
fields) \Q(u,u)\ < tj\ru\\ru\. This inequality in turn is deduced from an algebraic identity. For our 
problem such an algebraic identity does not exists for quadratic terms. However by analyzing the 
space time resonant frequencies and the null structure, all cancellations can be taken into account 
and this makes it possible for us to close our argument. Another difference is that here we need to 
use a cascade of energy estimates with controlled growth bounds. 



2. Sketch of the proof 

2.1. Ideas of the proof. 

1 ) Weighted Energy estimates. The basic estimate is the conservation of energy 

/ \v\ 2 dxdz+ [ (^l + \dh\ 2 -l)dx = [ ^G{h)i)dx+ [ ( y/l + \dh\ 2 - 1) dx 
J@ Jr 2 Jk 2 Jr 2 

which is a consequence of invariance under time translation. Higher energy estimates are derived by 
using the geometric structure of the problem, as was done in |29j. This essentially propagates the 
regularity of u at a rate of 3/2 spatial derivatives at a time, which can be seen from the model linear 

3 

problem dtu = —A?u. To simplify our presentation we elect to propagate the spatial regularity at 
twice that rate, i.e., 3. 

Since weighted energy estimates are derived from invariances of the equation |23j, we note 
that ransformations that keep (|WWj) invariant are: a) space translations, (h(t,x),ip(t,x)) i-> 
{h(t, x + 8), ip(t, x + 6)) for 5 £ M 2 ; b) space rotations, (h(t,x),tp(t,x)) i-> (h(t, Rgx),ip(t, Rgx)) 
for any 6, where Rg is the rotation of angle 6 around the origin; and c) scaling, (h(t,x),ip(t,x)) i-> 

^h(X 3 / 2 t,Xx),^(X 3 / 2 t,Xxf) for any A > 0. The vector fields associated with these transfor- 
mations are given by d, Q, and S. 

Since these vector fields are derived from the symmetries of the equation, they have good com- 
mutation properties. Combining them with the physical energy, this naturally leads to the weighted 
Sobolev spaces defined above. We show that the norm with the highest number of derivatives or vec- 
tor spaces grows slowly in time (like a small power of t): ||yl 1//2 uo|| H/ 9/2,2^ R2 ^ + H^l^Woll^o^^) ^ t s ' . 

These energy estimates are performed in Section [3j 



2) Dispersive bound. The weighted energy estimate relies crucially on the decay of A 1 ^ 2+a u in W K 
at the rate \jt. It is obtained from the linear estimate 

A^e-* tA3/2 f < jEIM^ 

L°° t II 
|fc|<3 

where Y(d) = \d\ L + for arbitrary small i, derived in Appendix [5j 



5,oo 



L 1 
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3) Space-time resonances. For the above linear estimate to yield the desired decay of A 1 l 2+a u in 
Wjl°° at the rate 1/t, it is easily seen that a uniform in time bound on A a u in W^ K suffices. This 
is obtained by the method of space-time resonances; the key analytic property of (WW) being the 
vanishing of the quadratic interactions at the space time resonant waves. We briefly explain this 
below. 

The space-time resonances method identifies the wave interactions which lead to significant 
contributions to the long time dynamics of solutions. It also presents an algorithmic method 
of how to deal with these interactions. It was introduced by the authors and proved very effi- 
cient in dealing with the global existence problem for a variety of dispersive equations: nonlinear 
Schrodinger [18] [19] . nonlinear Klein-Gordon [16], gravity water waves [20j . Euler-Maxwell |17j . 
We refer the reader to [15] for a comprehensive presentation of the method. 

For ()WWp the most significant interaction comes from the quadratic terms which can be com- 
puted explicitly by using the expansion for G{h) contained in Sulem and Sulem [30] 



d t h- 



Atp- 
■ Ah 



V • (/iV^) — A(hAip) + cubic and higher order terms 
^iV'i/'l 2 + ^|^1^| 2 + cubic and higher order terms. 

itA 3 / 2 



By writing the solution u = A 1 / 2 h + iip, in terms of its profile u(t) = e~ ltA "'" f{t), and using 
Fourier transform, one quickly discovers that the worst quadratic interaction, as far as space-time 
resonances are concerned, is u with u. The space time resonant frequencies in this case are given 

by = {(£,T])i £ = 0}. The vanishing of these quadratic interactions seems needed, which is 

the case for (jWWjh These calculat ions are carried out in Section [4l 

2.2. Detailed plan of the proof. Since local well-posedness is not an issue, the whole proof 
consists of a global a priori estimate. Define 

(1) The energy norm 



I energy 



2K 

d = supj^(t)' 

* 3=0 



(3+1)5 



A 1 ' 2 



II 



3(2A--j)+9/2,2 



i ,. def 
M| decay = SUp 
t 



2) The decay norm 

(t)~ s ' \\Y{dfA a u 

\t) l - & '-¥ 
1/2 L 

(3) The global norm 



+ ||r(<9) 2 yl a u| 



A+K 



+ sup sup 

t 0</3<l/2 



Y(dfA 



def 



vy 4,+K 



Y(d)A 



\^\\X — 1 1 U 1 1 energy 1 1 ^ 1 1 decay 

We can now define all the parameters appearing in the proof: 

• K > 10 controls the number of derivatives (see (|1.4I) ). 

• i « a,e: auxiliary parameter appearing in Y(^), used when summing dyadic pieces. 

• a* << 1: (fractional) number of derivatives of uq in L 2 for small frequencies (see (jl.4p ). 



a 



a* + 3l < < 1 : used in the definition of 



I energy 



and 1 1 It 1 1 decay 



• e << 1: controls the size of the initial data in X (see (jl.4p ). 

• R >> 1: large parameter which we do not set yet. It will be such that the solution satisfies 
\\u\\x < Re. 

• 5 = sup(Coi?e,a) << 1, for a big enough constant Co: time growth rate, appears in the 
definition of |M| en ergy 

• 5' = (2K + 1)5 « 1: time growth rate, appears in the definition of ||tx||decay 
The steps in proving the global a priori estimate are as follows: 
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Control of the energy norm. This is obtained in Section [3j There it is proved (Proposition 13.71 and 
Proposition 13.11")) that if the data is such that (|1.4[) holds and \\u\\x < Re, then 

IMIenergy - ^I 6 ' 

for some constant C[ uniform in R as long as Re « 1. 

Control of the decay norm. This is obtained in Section There it is proved (Proposition 14. 1| ) that 
if the data is such that (|1.4p holds and \\u\\x < Re, then 

Nldeca y <Cl(^+Nli)> 

for some constant C" uniform in R as long as Re << 1. 

Control of the \\ \\x norm. From the two previous points, we deduce that if the data is such that 
(|1.4p holds and \\u\\x < Re, then 

\\u\\ x < Ci (e+ ||u|||) , 

for some constant Ci > 1 uniformly in R as long as Re « 1. For e sufficiently small, the inequalities 

£ < Ci(e + x 2 ) and x > 

hold on two connected components, [0,xi] and [x2, oo] with x\ < X2 and xi ~ C^e, for a constant 
C<i- Picking R > C2, we get the desired estimate \\u\\x < Re by a continuity argument. 

2.3. Notations. We adopt the following standard notations 

• A < B if A < CB for some implicit constant C . The value of C may change from line to 
line. 

• A ~ B means that both A < B and B < A. 

• If / is a function over K then its Fourier transform, denoted /, or is given by El 

/(0 = ^/(f ) = ^ J e~**m dx thus f(x) = ^ J e^f(0 df. 

• The Fourier multiplier with symbol m(^) is defined by 

m(-d)f = F- 1 [mFf] . 
i 

• In particular, we denote (t being a sufficiently small real number) 

A d = f \d\ and d = A L + A~ L 

• The bilinear Fourier multiplier with symbol m is given by 

T m (f,g)(x) d = F' 1 j m(£,r])f(T])g(g - rf) drj. 

• The japanese bracket (•) stands for (x) = y/l + x 2 . 

^In the text, we systematically drop the constants such as \dii since they are irrelevant for our purposes. 
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3. Energy Estimates 
The conserved energy of system (jl.ip is given by 

(3.1) -^Diivsicai = / \v\ 2 dxdz + / 1 -_dS 

which is sufficient to give Hz(33) on the potential and H l {33) on h for small dh. Higher H s 
bounds and weighted norms bounds will be derived by commuting vector fields with the evolution 
equation of the mean curvature k of 33, and from Euler's equation evaluated on the boundary dU 
(see [29] for the derivation) 

(3.2a) D t n = -A@v ■ N - 2U ■ (V T v) 

(3.2b) N ■ D t v = -Nk + N ■ VAq 1 div(v ■ Vv) 

(3.2c) D t {z - h) = v 3 - D t h = 0, 

where D t = dt + v ■ V G Tr x t \dU is the material derivative on dU, A@ is the surface Laplacian, n 
is the second fundamental form of 33, N is the Dirichlet Neumann operator of Ql, T and _L are the 
tangential and normal component to 33 of the relevant quantity, Aq 1 is the inverse Laplacian on 
*3l with Dirichlet data. Here we also introduced the notation A\ • = trace (Ai^D), for two 
matrices, where A* is the adjoint of A. 

In order to compute commutators of vector fields with the above system we need to introduce 
some geometric notation. 

3.1. Geometric notations. Recall that U = {(t,x,z); (x,z) G denotes the space time fluid 
domain. Let 

dt = d t + (d t h)d z , Bi = di + {d l h)d z , i = i,2, 

denote vector fields defined on IA that are tangent to dlA. For any function <p defined on IA let ip^ 
denote its value on U, i.e., <pb(t,x) = tp(t, x, h(t, x)), consequently 

d a ip{t,x,z)\ z=h{t ,x) = d a (<p b (t,x)). 

Also note that since di G T33 are linearly independent, they can be used to construct an orthonormal 
frame on T33, which we denote by {e\,e2). 

For any vector e G T33 the covariant derivative T> e on T33 can be defined in terms of the outward 
unit normal iV _L T33 as follows. Writing V e = e • V for the directional derivative in M 3 , we have 

V e w = (V e w) T = V e w - (Vew) 1 - = V e w + (w ■ V e N)N, w G T33. 

For a vector- valued function on R 3 / = (/ J ), we also denote Df for the partial derivatives matrix 

Df = !/;,/'),.,. 

For any orthonormal frame {e a } of T33, the second fundamental form n, the mean curvature k of 
33, the Hessian T> 2 and the surface Laplacian A : a§ are given by 

n(e a ) = V ea iV, K = V ea N-e a , Vcj)(X, Y) = V x Vy<£ - V Vx Y<f>, 

where we sum over repeated indices. Note that using the frame {e a } the equation for k can be 
rewritten as 



(3.3) 



D t K = -Aggv ■ N - 2(V ea N) • (V ea v). 
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The harmonic extension to Q of a function / defined on SS is denoted by and the Dirichlet- 
Neumann operator by M 

W) = Viv/w : # -> M. 
We refer to Section [6] for the definitions and main properties of the spaces H S (S§) and H s (2t). 

3.2. Vector fields. Recall that we defined 12 and 5 on M 2 by 

(3.4) Q = x 1 d 2 - x 2 di = u% and S = + x 1 d 1 + x 2 d 2 = -td t + x%. 

For any function ip : U — > R with = </?|aw we have 

^tdt + x l <9^ <p 6 = (%tdt + + zd 2 + (Sh - h)d z \ <p 
A d 2 - x 2 di) tp b = (x 1 d 2 - x 2 di + (x 1 d 2 h - x 2 d 1 h)d z ) p 



m 
au 



Thus with a slight abuse of notation, (S, f2), can be considered as vector fields defined on dU by 
writing 

s = hd t + x % + zd z + zd z , z = (hd t + x%)h - h, 

Q = x l d 2 — x 2 d\ + (x 1 d 2 — x 2 di)hd z . 
The coordinates of the space parts of S and fl are denoted respectively (S l ) and 











\ s 3 ) 





x 1 > 










x 2 


and 




H 


1 S ) 


z + Sh-h ) 








V Oh I 



These vector fields can also be extended harmonically on 3! (S%, fly.), and thus on U: 

3 

S-h = 2 t9 t + x%a i + zd * + Z Hdz, % = x 1 d 2 - x 2 d x + lj^8 z . 

3.3. Commutators. Since on dU, S,f2,D t G TdZY, then [S,D t ], [fi,D t ], [D,S] G TdW, and 

(3.5a) [S,D t ]<p + ^D t <p= (^v* + Sv^Bnp 

(3.5b) [fl, D t ](p = Qv%<p - v x ~d 2 p + v 2 d x p. 

Remark 3.1. One can think of these commutators in the following way. 

• [S, Dt] + |-Df is a first order differential operator with coefficients depending on (v, Sv, dhv, dhSv) . 

• [f2, Dt] is a first order differential operator with coefficients depending on (v , Qv , dhv , dhQv) . 

To compute commutators of (D t ,S,Q) with the surface Laplacian Agg, we proceed as follows. 
Writing X = a(t)dt + X ■ V G TdU for any of these vector fields, we first note that for any vector 
field e G TS$ we have [X,e] G TSS since [X,e] G TdU and has no dt term. Thus 

[X, e ]=Xe- V e X = (Xe) T - (V e X) T 

We also note that for any orthonormal frame {e a } of TSS that depends smoothly on t, we have 

(3.6) (Xe a ) T -e b + e a - (Xe b ) T = X(e a ■ e b ) = 0. 
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With these observations we compute 

•VA /V = X (V ea V eo ^ - V Vea e a <p) 

= A 9 Xtp + 2V 2 (p(e a , [X, e a }) + V w ip 

= A<%X<p - 2V 2 <p(e a , (V £a X) T ) + V#<P, by §3M) 

where W = T>\ Xea \e a + V ea [X , e a ] — [X, V ea e a ] E TBB. By keeping track of the tangential terms we 
get 

W = {V [XM e a + V ea [X,e a ] - X(V £a e a - kN) + V Vea£a X) T 

= (V^,^ + V ea Xe a - V ea V ea X - X(V ea e a ) + X(kN) + V^ ea X) T 
= - (A m X) T + kXN. 

Thus 

(3.7) [Af, A^]^ = -2V 2 <p • V T X - (A@X) • (V T p) + «vfiV • (V T (^), 

and the commutators of Dt, S, and with Aa% are given by differential operators 

(3.8a) [D t , A<%]ip = -2V 2 ip ■ V T v - (A a v) ■ (V T \p) + nD t N ■ (V T <p) 

(3.8b) [S, A m ]y = -2A@ip - 2V 2 ip • V T (Zk) - (A&S) • {V T tp) + kSN • (V T v?) 

(3.8c) [Q, A<%]<p = -2V 2 ip ■ V t lo - (A m uj) • (V T p) + kQN ■ (V T <p) 

Remark 3.2. Using the relation Dfh = v% we get 

• [Df,Agg] is a second order operator with coefficients depending on w\ d = (d 2 h, (dh) 2 , Dv, D 2 v, 

dhDv,dhD 2 v). 

• [S, A i og\+2Aag is a second order operator with coefficients depending on W2 = (d 2 h, (dh) 2 ,Sdh,Sd 2 h). 

• [Q, Agg] is a second order operator with coefficients depending on w% d = (d 2 h, (dh) 2 , Qdh, Qd 2 h). 

To compute commutators of (Dt,S, Q) with the Dirichlet-Neumann operator M we proceed as 
follows. Again writing X = a(t)dt + X ■ V E TdU, for either Dt, S, or Q, we have for any g : U — > R 
and if : dU -> R 

(3.9a) #ttV S = VAT W5 - (£>X«)*V 5 

(3.9b) Xuvh = (Xip) H + A^2diy((DX H )*Vip H ) 

(3.9c) ^A ^ = A 1 (X n g) + A^2div((DX n )*V A^g) 

This gives the identities 

[X n ,VA^dw]g = (DXnYVA^dwg + 2VA ' 1 div ((£>X«)*VAq ^v^) + VA^ 1 (DX n • £> 5 ) 

and 

(3.10) [X,M]tp = (XN) -V m -N- ((DX n )*V m ) + N • VA 1 2div(( J DX w )*V^). 
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(3.11a) 
(3.11b) 

(3.11c) 



[D t ,M]ip = D t N ■ Vip n - N ■ ((Dv)*V<p H ) + N • VA l 2div((Dv)*Vip n ), 

[S,J\f]ip +Mip = (SN) -V m -N- (VZ n )d z tp n + N • VA 1 2div((VZ w )a z ^) 

= (ON) ■ V<pH - + N 2 d Wu - N ■ (V^)d zm 

+ N- VA - 1 2div((Va;| i )a 2 ^); 

The only term that prevents the order of these commutators to be apparent is N ■ VAq 1 2div. 
However by proposition (|6.3p 

N • VAg Miv : H s+ ^{9) _> H s (m) 
is a bounded operator. Thus if $ : H s (83) — > H s+ 2 (3>), is an extension operator, then 

N ■ VA^div^ : H s {£g) ^ H s {£g), 

is a bounded operator. 

Remark 3.3. Using the relation D t h = v% we get 

def 

• [Dt,N] is a first order operator with coefficients depending on w\ = (Dv,dhDv). 

def 

• [S, J\ \ + M is a first order operator with coefficients depending on W2 = (dh,Sdh). 

• [£2,J\f] is a first order operator with coefficients depending on d = (dh, Qdh). 
Bounds on commutators are given in the following lemmata. 

Lemma 3.4. Assume (v,dh) are smooth and \\(v, dh)\\ w ^ !00 < 1, for s > 2, then 



(3.12a) 
(3.12b) 
(3.12c) 



[D t ,V] ( p\\ L ,<\\(Dv,dhDv)\\ L »,\\V<p\\ L 3 

[D u V]<p\\ H . < \\(Dv,dhDv)\\ w ./ a , 00 \\V<p\\ H ' + \\Vcp\\ Ws/2 ,oo\\(Dv,dhDv)\\ H s 
[A,A/"]v?|| L 2 < ||u)i||l°°||V<^|| L 2, 

[D t ,N]<p\\H° < ||*l|liy[V2],oo||V^||^ + ||*l||^ ||V^|| VK [ s /2],oo. 

P,Ar]<p\\ L 2 <\\d 2 h\\ L ~\\v<p\\ L i, 

[V,M](p\\ B . < \\d 2 h\\ wls/2loc \\Vip\\ H s + \\d 2 h\\ H s ||V^|| W [V2],cx>. 



(3.12d) \([D t ,M]cp, <p)\ < [\\Dv\loo + \\Dh\\ w i,<*, WDvWwi^ + H^IU] \\<p\\ 2 H i/2 

Proof. [Dt,"D] is a first order differential operator whose coefficients depend on (Dv,dhDv), thus 
(|3.12ap is straightforward. 

To compute the norm of [Dt,M], given by (|3.11a|) . on H s we have to bound terms of the form 
G ■ V T , gAfip, \N ■ VA " 1 div(AV-), where G, g, and A are given in terms of W\ = (Dv,dhDv). By 
proposition 16.31 these terms can be bounded by 

||G-V T ^|| L2 + ||^|| L2 < \\(G,Gdh,g,gdh)\\ L ~\\V<p\\v 

\\G- V T <p\\ H s + \\gtf<p\\H> < \\(G,Gdh,g,gdh)\\ w[ s hoo \\V<p\\ H s + \\(G,Gdh,g,gdh)\\ H s\\V V \\ wmi0a 
\\N. VAoMiv^V^H^ < ||A|| L <»||V V || La 

\\N ■ VAq Miv^V^llH- < (1 + \\dh\\ wl s_ lhoo )\\AV^\\ HS + \\dh\\ HS -i \\AV<p\\ wm>oa , 
and this implies (|3.12b|h Equation (|3. 12c|) follows from a similar computation. 
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Finally to derive (|3.12d|) we note that ([D t ,N~\(p, ip) can be written as a sum of three terms which 
can be estimated as follows: 



ipG- V T ipdS 
gipMip dS 



<(||V T G|U» + \\G\\ L ~\\d 2 h\\ L ~)\v\ 2 L2 , 

A Jgs Jscd 

S3 

< [ll^lloo + ||-D/l|| W l,oo||Du|| W l j0 o + HD^IIco] M#l/2(^) 

where the last inequality sign follows from Proposition 17. 1\ and finally 
J LpN ■ V A l di\(DvVLp H )dS = j (p H dw(DvV(p H ) dx 



Jn-(d 



vV(p n )<p n dS 



+ 



< (\\Dv\\ W l,<x> + ||5 2 /l||Loo||Du||iao)||^||^. 1/2 , 



which implies the stated inequality 



□ 



Using these commutators we can derive a second order evolution equation for k. We start by 
choosing an orthonormal frame {e a } which is parallelly transported along the particle path i.e., 

D t e a = (V ea v ■ N)N. 

Applying D t to (|3.3|) . we obtain (see equation (3.15) in [2"5] ) 

D 2 k = -D t A m v ■ N - A m v ■ D t N - 2(D t V ea N) ■ (V £a v) - 2(V £a N) ■ (D t V ea v). 

Using the fact that 

D t N = -[(Dv)*N] T , and [D t , V e J = -V (VeaU) x, 
and the commutator formula [Dt,A&], we obtain (see equation (3.16) in |29j) 
D 2 k = - A m D t v ■ N - 2V ea N ■ V Ca D t v - 2A^v ■ D t N 

AV ea {D t N) + 2V {VeaV) rN)) ■ V ea v - K \D t N\ 2 - 2D t N ■ V {DtN) N 

Using Euler's equation evaluated on 88 to substitute for Dtv on 8$, we obtain 
D 2 k =A,u S N(k) - Vk h ■ A,%N - A^VAq 1 div(v ■ Vv) ■ N 

2A m v ■ D t N - (4V ea (D t N) + 2V (Veol)) T iV)) • V ' ea v - n\D t N\ 2 
2V ea N ■ (V ea VA^ l div(v ■ Vv)) - 2D t N ■ V (DtN) N 
' A^M(k) + R . 

Note that in terms of regularity on the boundary 88, we have from (|3.2p 

v ~ Dth, Aagv ~ D t n, N ■ D t v ~ Mk ~ D 3 h. 



(3.13) 



def, 
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Thus the nonlinear terms of Rq are 

Vk h ■ A^N = (d 3 h) 2 + cubic 

A^VAq 1 div(v ■ Vv) ■ N = DvD 2 v + cubic 

A@v ■ D t N = D 2 vd t dh + cubic 

V ea (D t N) ■ V £a v = d t d 2 hDv + cubic. 

def 

If we introduce w = (v,dh), then we can write 

(3.14) #0 ~ D t KD t dh + (8k) 2 + cubic ~ (Dw, D 2 w) 2 + (wf. 
Energy estimates will be derived using equation (|3.13p . 

Remark 3.5. The energy and weighted energy estimates can be done simultaneously. However we 
elected to separate the two for the benefit of the reader. The H s estimate is derived by commuting 
A = AogM with the equation which leads to a schematic equation d 2 Ag—d 3 Ag = {dAg)dg, while the 
weighted estimates which are derived from commuting S and I? lead to d 2 Sg — d 3 Sg = (d^g)dSg. 
Thus energy estimates for Sg require estimates on d s g. In terms of regularity, S acts therefore 
as 3/2 derivatives. To avoid fractional derivatives in energy (and weighted energy) estimates we 
estimate derivatives in multiples of 3. 

3.4. H s (33) estimates. For the remainder of this section all norms of v are computed for v\gg, 
thus H s stands for H s (3$), etc. Bounds of dh and v in H s will be derived from the equation 

D 2 n - A m Mn = R on dU. 

and as such they are implicit in [29]. 

Energy estimates are usually obtained in a straightforward manner by differentiating the equation 
and multiplying it by DtK. However in this case, there will be commutators present, and bounding 
them requires some care. 

Consider smooth solutions of 

(3.15) D 2 g - A /\'<i = F on dU. 

For functions defined on 3B let ( , ) denote the inner product on L 2 (3§), and define 

e g (t) = (VAfD t g, T>MD t g) + (AfA^Afg, A^Mg) = \\VAfD t g\\ 2 L2 + \\A^Afgf ' . 

Lemma 3.6. Smooth solutions to the above equation satisfy, if (Dv,dh) is small in W^ ,oc , 

e g < \\(Dv,dhDv,d 2 h)\\ w3 ,°o(e g + \\D t g\\ 2 L2 + \\Vg\\ 2 , ) + \\VNF\\ L 2^ 

ti 2 

Proof. Recall first that 

D t dS= (v ± n + V-v T )dS. 
Multiply ([3.15P by N A^N D t g and integrate over 38 to get 

e g (t) ={(v ± K + V-v T )VMD t g,VMD t g)+2{[D t ,V]MD t g,VMD t g)+2{V[D t ,M]D t g,VMD t g) 
+ {[D t ,M]Aatfg, A m Ng) + 2(M[D t , AaM]g, A m Ng) + (VMF, VND t g). 
Thus to estimate e g we use the commutator bounds given in Lemma 13.-41 These bounds imply 

\{[D t ,V]ND t g,VND t g)\ + \{V[D u N}D t g,VND t g)\ < \\(Dv,dhDv,d 2 h)\\ w ^\\D t g\\ 2 H2 , 
and by (|3.12dl) of the same lemma we have 

\([D t ,Af]A^g, AggNg)\ + |(AT[A, A^A% A m Ng)\ < \\{Dv,dhDv,d 2 h)\\ w3 , ac \\Vg\\ 2 H5/2 . 
Consequently 

e g < \\(Dv,dhDv,d 2 h)\\ wlt00 (e g + \\D t g\\ 2 L2 + WVgf^) + \\VAfF\\ L2y /%, 
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which implies the desired result. □ 

H s estimates for solutions of (|1.1|> follow in a straightforward manner from Lemma 13.61 Let 
n n = A n n, where A = Aogftf, then 

(3.16) D 2 t K n = AajM(K n ) + R n , 

where R n is defined inductively by 

R n d = f -[A, AVn-l +ARn-l = D t [D t ,A]K n -l + [A, A]D t K n -i + ARn-U 

with Rq given by f)3. 13|) . Thus 

n 

(3.17) Rn = J^-* (A[A,^]«i-i + [D t ,A]D tKi -x) + A n Ro 

i=i 

Define the energy 

£»(*) d - ^{VMD t K n ,VMD t K n ) + i(AA^K n ,^K n ) 
where ^4 = A<jgN. Finally, let 

n 

En(t) = / u £k(t) + -^physical) 
fc=0 

where n is a large integer and -^physical is given by (|3.ip . Then is equivalent to 

E n ~ ||AK||^3n+2 + 11^11^+7/2 + ll^Vllia + l|9/i||| 2 

~ II^Vll^+i + ll^ll^ +9/2 

Global energy estimates will be derived under the assumptions 
(E 2K {t) < e, for e < 1, for K large, 

\b(t) = \\d 2 h\\ w 3K+7/2,a a + 1 1 DV 11^3^+3,00 < tt~ l , for t>l. 
Which follow from 1 1 U 1 1 decay ^5 e - 

Proposition 3.7. Under assumption (|aip . solutions of (|3.13p u>ii/i data (ho,vo) 

E 2 k{0) ~ 1 1 «o 1 1^+4 + ll^o||^6if+9/2 < e 2 < 1, 

(3.18) £ 2 KW<e 2 i 2<5 , 
Proof. Recall that on the boundary 3S, we have from (|3.2|) 

v ~ D t h, A@v ~ A«) N ■ D t v ~ Nk ~ D^h. 
From Lemma 13.61 we have for n < 2i"T 



(al) 



(3.19) £ n < ||(^,5/ l I?i;,5 i / l )|| w i,oo(f n + ||A«n||£a + II^K i) + IpAARnllW^- 

From Lemma 13.41 and (JaT|) , we have for n > 4 

WVNA^DtiD^A}^-^ + \\VMA n - i [D t ,A]D t K i _ l \\ L 2 < - f ^K, 

where we bound the low derivatives in L°° and the high derivatives in L 2 . From (13. 14ft we know 
that Rq ~ DtuDtdh + (Ok) 2 , and thus 

\\VMA n R \\ L2 < jy/K, 
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by assumption ([aTj) , Consequently (|3.19|) can be rewritten as £ n < jE n , and by summing over 
n < 2K we get the stated inequality. □ 

Note that in terms of (ip,h), where ip is the boundary value of the potential ifa, Eq bounds 
%j) G H2(R 2 ) and h G H 1 ^), and thus Proposition [377] gives control on 
E 2K ~ jRr+9/2 + ~ \\A^uf 



|#6K+9/2( K 2) T ||l-"<-||jy6K+9/2( K 2) r ~ I K 1 ' « 1 1 H 6K +9/2 

where u = A 1 l 2 h + i^>. 

3.5. Weighted energy estimates. Weighted energy estimates will be derived by commuting 
(5, J?) with the k equation (|3.13p . The main difference, and difficulty, between the energy estimate 
and weighted estimates is the following. Commuting Dt with M gives a first order differential 
operator with coefficients depending on w\ = (Dv, dhDv), while commuting S or Q with J\f gives a 
first order differential operator with coefficients depending on W2 = (dh, Sdh), The crucial difference 
is that w\ decays ~ j in L°°, whereas w\ decays ~ ^7 (compare assumptions (jalj) and (|a2j) for 
instance). For this reason weighted estimates are slightly worse than regular energy estimates. We 
also would like to remind the reader that since S in commutators acts as 3/2 derivatives, we will 
treat (5, f2,A), where A = AggN, on equal footing; these operators are collectively denoted by r. 
In the sequel we will let D stand for either d or M, and X stand for 

X = S,Q. 

Finally, we define the spaces, 

(3.20) W k ' p (<%) = {ge W k ' p (<%); Mg G W k ~^ p {M)} 

(3.21) L*{&) 

(3.22) W k ' p (^) = {ip: 

i=0 

for p = 2 or 00, k,£ G N. For p = 2 we allow k G N/2. Note that W k > 2 (&) = W k ' 2 (^), however 
W k '°°{^) ^ W k '°°{m), due to the presence of M. 

The relation between v, h, and k in W k,p norms is as follows. 

Proposition 3.8. Under the assumption 

\\dh\\ w z ca + Hull^oo <C 1, for K > 1, 

solutions to system (|3.2|) satisfy 

(3.23) \\dh\\ w a+i,2 [aJ) ~ INI^^ + \\dh\\ L 2 K{m 

(3-24) 1^11^+2,2^ ~ \\D t K\\ w a^ m + M\ w «g {m) + \ \ dh 1 1 L 2 k {a) + || v L \\lI k {3S) ■ 




Proof. The first inequality follows from k = — (1 + |5/i| 2 ) 2 A@h, commuting with A#, and 
standard elliptic theory. To prove the second inequality we note that since v is the gradient of a 
harmonic function, then by equation (16.3P 

\\ D2 v\\ L 2^ < \\v\\ H 2^ < \\v < \\Ad&V \\ L 2(&) + ||i'" L ||l2(^) 

From equation (I3.2ap 

Aggv 1 - = -D t K - v^lUl 2 + (V ■ U)(v T ), 
the commutation relations of X , and by the fact that | XT | + |X?n| < \dh\ + \D 2 h\ + \D 3 h\, we conclude 
the stated result. □ 
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Set hj = K2K-j = A 2K ~^k, for any < j < IK. It solves the equation 
(3.25) D^kj = A<%N(kj) + Rj, 

where Rj = RiK-j is defined by (|3.17p . Next, let 

where I = . . . 2K — j. Finally, define the weighted energy that will be controlled: 

n 2K-j 



Fn — ^ ^ ^ * ^ j ^physical* 



j=0 e=o 



Notice that 



F n (t) ~ \\D t K\\ 2 w3(2K . n)+2 , 2 + ||«|| 2 3(2K _ n)+ 7 i2 + U 1 / 2 ^ 2 + \\dh\\ 2 2 



\v\\ 2 w 3, { 2K-n)+A,2 + \\dkf 3(2jf _ n)+ 9 >2 + H^Vlli 



2 

3(2A'-n)+§,2 • 



To derive the weighted energy estimate we need to commute X with the operator D 2 — AggN- 
The next lemma illustrates how and where terms such as Sdh appear in the commutators. 

Lemma 3.9. Let g denote a smooth function on dlA and write 
(3.26) D 2 t g-A m Mg = F on dU. 

Then 

(3 27) WVmX^xN- D 2 ]g\\ L2 < (||^||^ 4i<JO + ||Aft;||^ iCO ) (\\Dg\\ H * + \\DD t g\\ H z) 

+ \\VMF\\ L 2. 

Proof. For X = S, we have from the commutation relations (|3.8f) and (|3.1ip 
[S, A m N - D 2 ]g = (A&[S,N] + [S, A<%]N - D t [S, A] - [S, A] A) g 

=A m {VsnOh ~ {MZ)d z g H + N ■ VA l 2div((V Z n )d z g n )) 

- (lA m Mg + 2V 2 Mg ■ V T (Zk) + (A<gS) • (V T A"c/) + kSN • (V T A" fl )) 

- A(^ + Sv l )d i9 - {\ l + Sv^diDtg + 3D t F. 

Writing A\jz H {d z g) = Ao$ (N ■ VAq l 2d\v{(y Z-^)d z gy^ (which is an operator of order 2), this 
commutator can be written as 

[S, A m hl - D\\g =A 3 {SDh)D 3 g + A 2 (SD 2 h)D 2 g + A 1 (SD 3 h)Dg 

+ Avz H (d z g H ) + B 1 (v,Sv)DD t g + 3D t F, 

where we exhibited the dependence of the coefficients on the least regular terms and terms with 
the least number of D. Here we introduced the notation A3D 3 = Ylf=i A\N l D 3 ~ l , and so on. This 
implies 

\\VN[S, A m M- D 2 ]g\\ L2 < (\\SDh\\^ tOB + ||St;|lw*,<») (\\Dg\\ m + \\DD t g\\ m ) + \\VMD t F\\ L ,. 
Similar estimates hold for X = Q. □ 
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Remark 3.10. Note that this lemma, when combined with Lemma \3.(A gives an estimate with 1/2 
regular derivatives loss on g. However since the loss is in regular derivatives, we can handle this 
by regular energy estimates. Moreover since terms like SDh decay at a rate less than t , this will 
cause a cascading effect on the growth rate of the weighted energy estimates. 

Bounds for F n (t) will be derived under the assumptions 

(F (t) < e 2 t 2S , \\v\\l 4 , 2 + \\dh\\* < e, for e « 1 



(a2) 



a(t) = f ||0/i|| W 9/ a ,« + IMI^oo < et^+ s , t > 1, 
b(t) = f Wd^W-y^ + ||^||^3 l0o < et-\ t > 1, 



I decay 

Proposition 3.11. Under assumption (|a2p . solutions of (|3.13p with initial data 

\\v \\ 2 w 4,2 + \\dh \\l r9/ 2,2 < e 2 « 1, 



>K 



satisfy 

(3.28) Fj(t) < e 2 t 2 (i+i)5 5 
for 0<j< 2K. 

Proof. Step 1: the ODE controlling Fj. The proof will be constructed inductively on j. The case 
j = is the non- weighted estimates. Since the difficulty in controlling Fj is in controlling high 
order weighted derivatives, we will only keep track of the highest order weighted derivative terms 
in Fj, that is J^j 1 . Recall that kj = A 2K ~^k satisfies (|3.25p . and that 

(3.29) &f(t) = ^\\VND t X j kj\\ 2 L2 + ^ A'A/A'.T'/V,. A^A'.T'/,-, . 

For 1 < j < 2K, X 3 kj satisfies the equation 

(3.30) D 2 X j kj = A^M(X j kj) + Kj, 

j 

(3.31) Kj d = f ^ X'-< ([X, A m M - D^X^kj) + XiRj, 

i=i 

and consequently from Lemma 13.61 we have 
(3.32) 
d 



j t ^j{t) =((v x k + V- v 1 )VND t kj,VMD t kj) + 2{[D t ,T>}ND t k j ,VND t k j ) + 2 (V [D t , A] D t kj , VND t 
+ (\D t ,N]A m Nkj, A :M Mkj) + 2{M[D t , A%N)kj, A m Nkj) + (VNKj ,T>N D t kj) , 

and 

< \\(Dv,dhDv,d 2 h)\\ W3 ^ + WVMllxWvyffi < ~&j + 
Thus to prove the proposition we have to bound Kj, which will be done inductively. 
Step 2: the case j = 1. For j = 1 

Ki = [X, A S N - D 2 ]k x + XRx 

has to be bounded carefully due to the presence of Xdh in the commutators whenever X = S or 
17. To estimate this term when X = S, we use Lemma 13.91 to obtain 

\\VN[S,A^N- A 2 ]^ilb < (WSDhW^ + \\Sv\\^\ (WDk^H, + \\DD t kx\\m) + \\VND t R x \\ L 2 
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From the regular energy estimates and assumption (|a2|) we can bound the right-hand side as follows 



\\Dki\\ Hi + \\DD t kx\\ H 2 < ^/E^k < et 5 
\\VMD t R 1 \\ v ><ly/E^< 



t V ~ £1-5 > 



Thus ||XW[«S, A.^TV- D t 2 ]Ki|| L 2 < p^-. Similarly, observe that 

2K-1 

Rx = RiK-i = A 2 *- 1 " 1 (Dt [Dt, A)Ki-i + [D t ,A]D t Ki^) +A 2K ~ 1 R , 

i=l 

is a sum of products of the following type: a function which carries the highest number of derivative, 
which can be estimated in L 2 ; and a function that decays ~ \ in L°° due to Assumption (ja2[) . 
Therefore, we get 



Similar estimates holds when X = Q. Thus we conclude 
which proves the proposition for j = 1. 



Step 3: the case j > 2. For j > 2 we proceed inductively in j to estimate ||£>.A/"7£j||£2. Assume 
we have verified the proposition for j — 1. To show it holds for j we have to apply X 3 to the kj 
equation. Recall that 

[X,A m M - D^X^Hj =A 3 (XDh)D 3 X i ~ 1 k j + A 2 (XD 2 h)D 2 X i ~ 1 k j + A 1 (XD 3 h)DX i ~ 1 k j 

+ Ayz^d.X^kjn) + B 1 {Xv)DD t X i - 1 k j + SDtX^kj. 
In equation (|3.3ip . these terms generate quadratic terms of the form 
A 3 (X e+1 Dh)D 3 X^ 1 k J +A 2 (X e+1 D 2 h)D 2 X 3 ^ 1 k j 

+A 1 {X i+1 D 3 h)DX j ~ i ~ 1 k j + B 1 {X i+1 v)DD t X j ~ e ~ 1 k j , 
(for < £ < j — 1) as well as terms of the type 

A, aJ (N • V^ l 2div({X^JX H )d z X^~ l k jH ) , < I < j - 1 



(3.33) 

A<g ( N • VAq 1 2^ 1 ( (VX W ) • {VXudzXi-^kju) ) ) , 1 < £ < j - 1, 



and finally terms which more regular and of higher degree. 

The first term in (|3.33|) is easy to bound, while the second term term gives rise to smooth regular 
cubic terms and to quartic terms of order —1 of the form 

(3.34) (a^n) • vAq x x l n l (jyx n ) • {V X H d z X^~ l k m ] 

Thus all these commutations lead to quadratic terms 

i-i 

Qi = s rA 3 (X e+1 Dh)D 3 X j - e - 1 k i +A 2 (X e+1 D 2 h)D 2 X j - e - 1 k i 

(3.35) j^o 

+A 1 (X i+1 D 3 h)DX j - £ - 1 k j + B 1 (X i+1 v)DD t X j - i - 1 k j 
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quartic terms of order —1 of the form 

(3.36) Kj = (A@ N) • VAq 1 2Xfr 1 ((VX M ) • (VX^*^- 1 ^)) , l<£<j-l, 

and more regular degree 3 or higher terms. Thus 

IZj = Qj + /Cj + X 3 Rq + similar or more regular terms of degree 3 or higher. 
To bound ||2?W72.y||^ 2 , we deal first with the quadratic terms Qj. In this case we bound each term 
in L 2 or L°°, depending on whether £>Kot£<K — 1. To estimate A^(X £ SDh)D 3 X : >~ £ ~ 1 kj , 
with I > K, we proceed as follows: 



\\VNA 3 (X e SD h)D 3 X j 



-t-lz. 



Kj\\ L 2 



< II D A X 3 



3 \\w 2 '- 



\X £ SDh\\ H 2. 



The above expression is bounded by et 1 \/F~j . 
For £ < K - 1, 

e 



\X e SDh\U 



< 



by assumption (|a2 



l W 2, oo f.\-8 ■ 

3 X j - e - 1 K j \\ H 2 < y/Fj-t-i < ei°'"^ )5 , by the inductive step. 



Therefore 



\\VNA^SDh)D^X^- l kj\\ L 2 < WX'SDhW^WD 3 ^- 1 ^ < -^js- 

The remaining terms in Qj are estimated in the same manner. The cubic terms are much easier 
to bound since they appear from multiple commutations of X and thus they have fewer derivatives 
and better t decay. 

The quartic terms ||XW7Cj||£2 can be estimated in the following manner: when T>j\f hits A^N 
we bound the term by 



(VMA m N) • VAq 1 ^" 1 l(VZ n ) ■ (VZndzX^^kjH 



L 2 



< \\dh\\ H , vAq ((vz H ) ■ {yz H d z x^ l k ]H: 



< \\dh\\ H A x^ 1 ((vz n ) ■ (yZndzXi-^Rju) 



By distributing the X^ 1 vector fields on the cubic term we end up with a highest derivative of 
Xi~ 2 , and at least one term with less than K derivatives. Thus 

CPA 'A /.V ; • VA^Xfr 1 ((VZ n ) ■ (V Z H d z X^-% H )) \^ < 

The remaining terms are quartic with fewer derivatives and are of positive order, thus their estimates 
are straightforward. 

This leaves \\VAfX J Rq\\ L 2 to estimate. Recall that Rq = R2K-3 is given by (|3.17p 

2K-j 

R2K-3 = ^ (D t [D t ,A]Ki-x + [D t ,A]DtKi-!) +A 2K ~ j Ro, 

i=l 

so these terms can be bounded by 

\\VNXiR4 L 2 < ~^Fj. 



Putting all of these bounds together implies 

d 



which gives the desired estimate after summing over j. 



□ 
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4. Decay estimates 

4.1. Rewriting the equation. Recall that the equations for the trace of the potential and the 
graph of the surface are given by (see (|1.2p ) 

d t h = Atp-V- (hVip) - A(hAip) + Ri 
d t i> = Ah- i|VVf + l\A^\ 2 + R2 

where Ri and R2 are terms of degree 3 and higher, and where A is the Fourier multiplier of symbol 
|£|. By introducing the variable H = A l l 2 h, the above equations can be written as 



(4.1) 



8 t H = A 3 / 2 i; - A l ' 2 V ■ (A-^HVip) - A 3 / 2 (A- 1 / 2 H Aip) + R x 
dtj, = -A 3 / 2 H - + \\Ail>\ 2 + R 2 , 



where R± = A l l 2 R\. The above system can be written in a compact form in terms of the complex 
variable u = H + tip: 

(4.2) d t u = -iA 3 ' 2 u + T m++ («, u) + T m __ (u, u) + T m+ _ (u, u) + R, 

where m±± are linear combinations of 

(4-3) def 

and mi((,(- 77), 

def 

and we set R = R.\ + iR 2 - The next step is to define 

and to write Duhamel's formula for /: 
(4.4) 

/(*,£) = £<>(£)+ E / /e^^^Wm^^^c^^^^e-r?)^^ 

(Tl,r 2 ) = (++), (+-),( ) ° 

+ f e is ^ i/2 R(s,0ds, 
Jo 



where /+ = f,f— = /, and the phases are given by 

<P±,±^v) = \a\ 3/2 ±\v\ 3/2 ±\z-v\ 3/2 - 

Sometimes, the exact structure of (|4.2p will not matter and we will write it in a simplified form 

(4.5) d t u = -iA 3/2 u + T m (u,u) + R, 

with m being a linear combination of mi and ni2. 

We will also often write indistinctively / for / and /; similarly, we will not distinguish between 
g is|D| 3 / 2 anc j e -is|D| 3 / 2 _ T/hjg alleviates somewhat the notations, and has of course no impact on the 
estimates. 

4.2. Examination of the bilinear symbols and the phases. 
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4.2.1. The symbols. Start with the two multilinear symbols 

rh^v) = f Z-(Z-v)-\m-ri\ 

def 

™2{£,v) = 7] • (£ - 77) + |r/||£ - 77 1 

Due to the relation 7711(77, £) = — rfi2{— — 77), it is equivalent to prove estimates for either of these 
symbols. 

Using the simple fact that 



\x\ 2 \\x\j V \ x \ J V ix 

we can deduce that for \rj\ « |£|, 

Mt, r?) = iei 2 - 7? ■ e - iei (lei + ^ + ° (jf 

= 0(|7 ? | 2 ). 

Therefore 

• If If I « M, then |7Tti(£,77)| < |C||r/| 

• If << If I, then |mi(£,r/)| < |r/| 2 

• If k~ »?l << M, then |7fii(£,r/)| < |f - 7/||7/| 
and thus 

. If |£|«H, then ^2(^77)1 <|f| 2 

• If |r/| « |f|, then \m 2 (Z,r])\ < \rj\\£\ 

• If If - ??| << hi, then |77i 2 (£,r/)| < |f - 7?||t?|. 

In terms of the classes defined in Appendix l8.2l this means that rh\{^,rj) G .A/f 2 ' 1 ' 2 ' 1 and 7712(^,77) E 

^2,2,1.1. 

Now 

|£|l/2 

^ = j^jT72 ^ aIld m2 ^' ^ = m2 ^' 

This implies that mi (£,77) £ .A/f 2 ' 3 / 2,3 / 2 ' 1 , mi (£ ; £ — 77) £ ^2,3/2,1,3/2 ^ and m2 (£ ;7? ) <= M 2 ' 2 ' 1,1 . 
Finally, since any m±± is a linear combination of mi (£,77), 77i2(f, 77), and mi(^ — 77), we deduce 
that 

m±± € A* 2 ' 3 ' 2 ' 1 ' 1 . 

4.2.2. T/ie phases. Recall that for a phase 0, the space, time, and space-time resonant sets are 
given by 

T = {(/> = 0} , S = {d r ,(/) = 0} and TZ = Sr\T. 
Let us study these sets for each of the phases <j>±± . 

• First 0++: it is obvious that 7++ = {(0,0)}, which is the most favourable case. We will 
therefore deal with this interaction by using a normal form transform, producing the symbol 
5±±, which belongs to ^1/2.8/2,1,1. 

• For 4>-\ , we have = +| ^172 + f |g-^/2 and hence we find 5^ = lZ-\ = {f = 0}. 

Thus in this case estimates will be obtained by an integration by parts in 77. Expanding 
d v (f)-\ around f = gives: 

^=^H( { -ra)ra) + °(i$' 

When integrating by parts in 77, the symbol will be changed to r^ d^tp, which, by the 

above estimate, belongs to ./Vt 3 / 2 ' 1 / 2 ' 1 ' 1 . 
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• Finally, for 
and 1Z = 



we have d n 



3 v 



•2^p72 + 2^p7 



3 



arid hence we find that = {£ = 2r/} 



{(0,0)}. As we will see, by splitting the (£,77) space into regions where and 

4_ and 4>a 



d v 4> do not vanish, this case can be reduced to estimates similar to those for $ 

4.3. Invariances and commutators. The solutions of equation (14. 2D are invariant by translation 
u i—)- u(x + 5), 5 £ M 2 , rotation u 1— > u{Rqx), and dilation u h-> -TT/2u(X 3 ^ 2 t, Xx). 
The generators of these transformations give the vector fieldsj 

^ = (di,d 2 ) 3 , Q = x 1 d 2 - x 2 di = u%, and S = -td t + x%. 
We shall also use the space part of 5, namely 

E = x%. 

The vector fields d 3 , f2,S are collectively denoted by r, and the multiindex notation for powers of 
r is used. 

Of course, d 3 commutes exacly with the linear part of the equation, and the linear group. As for 
the nonlinear part of the equation, it can be expanded into pseudo-product operators (by translation 
invariance) and regular derivatives can be commuted with pseudo-product operators according to 
the Leibniz rule. 

The commutation of the vector fields i? and S with the linear part of the equation and the linear 
group are given by 

n, d t + iA 3/2 \ = 



(4.6) 



S, d t + iA 3 ' 2 
f2,e ±itA3/2 



5, e 



±UA 3 / 2 



iA 3 ' 2 



0. 



Next, to commute these vector fields with the nonlinear part of the equation, expand it into a 
series of multilinear operators: in other words, write (|4.2[) as 



d t u 



iA 3/2 u + ^2L k (u,...,u), 

k=2 

where is /c-linear. The invariance by translation and rotation gives immediately that the Leibniz 
rule applies to d and ft, namely 

8L k (u, ...,u) = L k (du, . . . , it) H h L k (u, ...,du) 

QL k (u, ...,u) = L k (flu, ...,u)-\ h L k (u, flu). 

As for the dilation operators, notice that the scaling invariance gives 



(4.7) 



Ll 



lj^u(X 3 / 2 t, Xx), ... , -±^ U {X 3 '\ Xx)) = X [L k (u, ...,«)] (X 3 /\ Xx) 



A 1 / 2 



Taking the derivative in A yields the modified Leibniz rule: 

SL k {u, ...,«) = L k (Su, ... ,u) -\ h L k (u, ... , Su) 



k 



+ 1 L k (u,.. .,u). 



Finally, we will need to commute vector fields with Fourier multipliers of the type A a . The formulas 
are easily computed: 

[d, A a ] = [a, A a ] = 
[S,A a ] = [S,yl a ] = aA a . 



Notice that the vector field S equals the generator of the scaling transformation up to an additive constant only. 
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In particular, we shall remember that 

F A a = A°T + { lower order terms } 
T,A a = A a T, + { lower order terms } . 



4.4. The estimates. We prove the following a priori estimate (recall that / = e li ^ 12 u and that 
a = a* + 3i). 

Proposition 4.1. Assuming \\u\\x < oo and (T77 



(4.8) 
(4.9) 

(4.10) 

(4.11) 
(4.12) 

(4.13) 



Y(dfA a u\\ wS , 2 <(e+\\u\&)(ty 



\\Y(dfA a Zf\ 



< (e+\\n\\ 2 x ){tf 



7+K 



Y ( d) 3 A l/2+a-0 



II 



< 

W 4+K 



£ + U 



X 



(t) 



1-S' 



\\Y(d) 2 XA a f 
Y{d)A l ' 2+a -Pu 



if* 

7,2 <e+IMIx 



t/0 < P < 1/2 



|y(a) 2 ^|| W 7, 2 <e+\\u\\ x 



< 

^5,00 - 



e + kt 



A' 



*/0 < /S < 1/2. 



This proposition gives in particular the desired a priori estimate 

IMIdecay ^ e + IMIx- 

Remark 4.2. Notice that one of the reasons why many derivatives are lost in this argument is 
that we cannot interpolate easily the derivatives S and Q, namely we do not have a Trudinger type 
inequality for such derivatives. 

4.5. Proof of (14. 8|h The growth estimate (14. 8p follows from the result of Section [3] as far as high 
frequencies are concerned. Thus it suffices to control Y(d) 3 A a T k u in L 2 for \k\ < K + 8. We fix 
such a k and apply Y(dfA a r k to (j4"3]) . This gives 

Y(dfA a r k (d t + iyl 3/2 ) u = Y(d) 3 A a r k T m (u, u) + Y(d) 3 A a r k R. 



(4.14) 

Let us discard for the moment the remainder terms; we will come back to them in Section [4.111 

Next apply Leibniz rule (see Subsection (|4.3p ) to commute r k with the linear equation and the 
bilinear term. Out of the many terms coming out, we single out two representative examples: 
Y(dfA a T m (r k u,u) and Y(d) 3 A a T m (r k / 2 u, r k / 2 u) (assuming for simplicity that k is even). Thus 
the equation now reads 

(d t + iA 3 ' 2 ^ r k Y(d) 3 A a u =Y{d) 3 A a T m {r k u, u) + Y{8) 3 A a T m (r k/2 u, r k/2 u) 

+ { similar terms } . 

To estimate the first term on the right-hand side, recall that m £ _A/f 2 ,3/2,i,i anc j uge Corollary 18.11 
to get 



Y(d) 3 A a T m (r k u,u) 



< 



< 



< 



r k A a u 



A 3 / 4 



U'2 



r k A a u 




A 3 /\ 


+ 


r k A a u 




A 3 I\ 






2 








ir- 




W 2 '°° 



\u\\x 

(t) 



r k A a u 



i /4.\S'—1 II ||2 

+ (*> \m\x- 
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To estimate the second term on the right-hand side, we also rely on Corollary 18. II to obtain 



Y{dfA a T m {r k ^u,r k ' 2 u) 



< 



< 



r k / 2 A a u 

u \\x 



t 



H 2 

r k A a u 



A 3/4 pk/2 



U 



U -2 



Performing an energy estimate, (|4.14p gives thus the differential inequality 
d 



\Y{d) 6 A a u 



(t) 



\Y(dyA a u\\ w o, 2 _ + (t) 



u 



x • 



The hypothesis on the initial data (|1 .4[) and Gronwall's inequality give the desired conclusion. 



4.6. Proof of (14. 9]) . It follows from (|4.8p in a very similar, but slightly simpler, way to how (|4.12p 
follows from (|4.1ip . We therefore refer the reader to the proof of (|4.12p . 



4.7. Proof of (HTTP]) . It follows from and Proposition [5J] 

4.8. Proof of (liTTT]) . 

4.8.1. First reduction. The estimate (|4.1ip is the crucial one, and the one for which space-time 
resonances will play a role. We want to control \\Y 2 (d)r k A a d^u\\2 for \k\ < K + 4 and |j| < 7, 
which are from now on fixed. We will denote 

z = r k d j u and g(t) = e itA3/2 z{t). 

Start by applying Y 2 (d)A a r k &' to the equation (|4.2p . Using Leibniz rule (see Subsection I4.3[) . 
r k di can be commuted with the linear and quadratic terms, which produces many terms. We only 
record the extremal ones: these are the most difficult to estimate, and we will focus on them. Thus 
we find 



d t Y 2 (d)A a z = - iA^ 2 Y 2 (d)A a z + Y 2 (d)A a T m++ (z, u) + Y 2 (d)A a T m __ ( 



z,u 



(4.15) 



+ Y 2 (d)A a T m+ _ (z, u) + Y 2 (d)A a T m+ _ (u, z) 
+ { Mixed terms } + Y 2 (d)A a r k d j R. 
Proceeding as in Subsection 14. 1\ this can be translated in Fourier space into 

Y 2 m\ a d(t,o 



TL,2=± 

{ Mixed terms } + [ e is ^ 3/2 Y 2 (£)\£\ a r k di R(s, ds. 
Jo 



The constant term ^ 2 (£)l£| Q ^o(£)> t>y assumption, can be bounded in L 2 by e. 

In the three following subsections, we analyze separately the terms corresponding to ++, 
and +— in the above sum to show that 



s|D| 3 / 2 







Y(d) 2 A a r k d j T m±,± (u, u) ds 



< 



u 



x ■ 



The last term, which involves the remainder R, will be dealt with in Section \4. 1 II where we show 



e ls W 3/2 Y(d) 2 A a r k d j Rds 



o 



< \\ufx- 
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4.8.2. The ++ case. Let us consider for now the term 



o 



y\m\ a / / e fa *++fr")m(£, T])f(s, V )g(s, £- V )dri ds 



where m € J\4. 2,i / 2,1,1 (recall that we denoted indistinctively / for / or / since the distinction is 

r-t 



irrelevant for our purposes). Using the identity e iS * = jzd s e tsr ^ to integrate by parts gives 



(4.i6a) Y 2 mr 

(4.16b) 
(4.16c) 
(4.16d) 
(4.16e) 



e <»*++K.i) m (f , n)f{8, rj)g(s, £-r])drj ds 



Y z m\ a e 



Y 2 (m\ c 



f(t,v)g(t,£, ~ »?) drj ds 



/(0,r?)<7(0,£- V) d V ds 



Y 2 (m\ c 



s <t> ++ &v) d s f(s,r,)g(s,Z- V )dr l ds 

^^vTTl^' W'ft 3 ' t-^dri ds - 



Keeping in mind that belongs to .A/f 1 / 2,3 / 2 ' 1 . 1 ; applying Corollary 18.11 gives 



|(1™1>|| 2 < 



^1/4 



u 1 



I A' 



The term (|4.16cp is easier to treat, thus we skip it. Proceeding as in the proof of (|4.12p below, we 
can prove that 



YidfA^dtf 



2 - t 

This gives, with the help of Corollary 18.11 



<™ and 



< \Mx 
1-6 



2~ <t> 



idHMiii. 



< 



sy l3/2 



syi 3/2. 



r/.s 



< 



JO Hi 



ds 



w 1 



\\u\\x \\u\\x 
<*) (s) 1 - 6 ' 



ds 



< \W\\x- 



Finally, the term (|4.16e|) can be treated in a very similar way, thus we omit the details here. 



4.8.3. The +— case. Let us consider here the term 

pi 



Y 2 m\ a / eW + -&v) m far,)f(8,ri)g(8,Z-Ti)dT,d8 
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where m G TV! 2 ' 3 / 2 ' 1 ' 1 . The piece of the integral corresponding to s € (0,1) is easily dealt with; 
thus we shall consider in the following 

y\i)\i\ a f I e is *+-^m(Z, V )f(s, V )g(s,{;- V )d V ds 



(the aim of this manipulation is to allow for the upcoming integration by parts, which will produce 



a - factor, not integrable near s = 0). Using the identity e ls 
parts gives 

(4.17a) Y 2 mr 

(4.17b) 



d n 4> ■ d ri e ls< *' to integrate by 



•"m(£, r)) f(s, Tj)g(s, i-rfidr) ds 



(4.17c) 
(4.17d) 



Y 2 m\ c 



Y 2 (m\ c 



is\d v <i) + _{i,ri)\ 2 



f(s,v)d{ s ,£- V) drj ds. 



Notice that, due to the vanishing properties of m and d v 4> established in Subsection 14.21 the symbol 
i\dZf-ltv)\ +1 \l-v\ belon S s to the class 7W 1 / 2 ' 1 / 2 ' 1 ' . Therefore, by Corollary EH 



10323)11 



(4.18) 



< 



< 



2 — 

1 1 

1 8 



Y 2 m\ c 



is\d v (j) + -(€,r))\ 



;f{s,v)d v g{s,C- i)drjds 



ir 1 



i|8^^+_(€,'7)t 2 le-'7l 



Y(d)e isA3/2 Axg 



w 1 



e isA3/2 f(s),e isA3/2 Axg(s) 



ds 



ds 



On the one hand, we can interpolate between I? and L°° to obtain 



/1A < ^\\A-u\\ w ^\\A^u\\ wha 



< 



mix 



( 4 - 19 ) ^1,4 ~ V " """" ~ y/(t) 

On the other hand, the pointwise bound 

leil^OI Z \d(0\ + \n g (®\ + \fg(o\ < + 1 + 1^(0 1 + * 

entails the L 2 bound 

\Y{d)A 1 l 2 Axg\ Hi < \\Y(d)A 1 / 2 g\\ 2 + \\Y(d)A 1 / 2 g\\ w ^ + t\\Y(d)A 1 / 2 d t g\\ 2 < (\\u\\ x + \\u\\ 2 x )(tf 

(the term involving the time derivative being bounded as in the proof of (|4.12j) below) and, by 
Sobolev embedding, 
(4.20) 



ir 1 



Y(d)e isAZ ' 2 Axg < Y(d)e isAi ' 2 A 1 ' 2 Axg = Y(d)A x l 2 Axg 



u 1 



i-i 1 



< 

f->*> 



W\\x + \\u\\ 2 x )(tf 



Plugging flUgD and P~20|) in the last line of (^18|) gives 

11(0373)112 < 



1 1 IMU n n / \<5' j <r n it 2 
u\\x{s) ds<\\u\\ x . 



S \/S 
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The term (|4.17c|) is bounded in a similar fashion. Finally, in order to deal with (|4.17d|) . observe 



that d v ■ 



(4.21) 

\\Y 2 (m\ a mm> 



belongs to ^V4 1 / 2 > 1 / 2 >°>° ) thus we can use Corollary 18.11 to get 



< 



< 



e isA3/2 T 



F 2 (£)ISI a <V 



sA 3 ' 2 



\\Y{d)z\ 



w 1 



A 1 ' 2 



i4 + \\Y(d)u\\ wlA 



u 



ir 1 



is\d v fa 


-(£,' 


?) |2 


, \ is/I 3 / 2 / 






9\ 


.)) 


AV 2 z 




ds 









f(s,ri)g(s,€- rj)dr]ds 
ds 



By Sobolev embedding, and using the bound for H-A 1 / 2 ?^^^ in (|4.19p . as well as a similar bound 
for ||yl 1 ' /2 2;|| H /i,4, the above can be estimated by 



\Y 2 m\ a mm 



< 

2 ~ 



< 



1 1 



S L 
*HI« 



Y{d)A 1 / 2 g 



H 1 



A l '\ 



ii' 1 



+ 



Y{d)A 1 ' 2 f 



H 1 



A 1 / 2 ; 



w 1 



ds 



x 



s \/S 



u\\x(s) ds<\\u\\ 2 x . 



4.8.4. The case. Upon partitioning the (£, rj) space into regions where cj>, respectively d^cj) do 

not vanish, one can proceed as in the ++ case (if 4> 7^ 0) or case (if d v (j) / 0). We refer 

to [19] for a similar partitioning. Indeed, recall that for £ = 2rj, <j)(2rj,rj) = (2 3 / 2 — 2)|r/| 3 / 2 and 
hence = {(0,0)}. We can partition the plan £,77 into two regions R 4 = £2\ U O2 such that 

I ^ 1 3/2 _|_ 1 ^ 1 3/2 < |0 J j n region i?i and |£| 1//2 + |??| 1//2 < \d v (f) | in the region 4?2. 

More precisely, we can take a smooth cut-off p : R + — > M + such that p(a) = 1 for a < 1 and 
p(a) = for a > 2 and then define 

(4.22) f&^m^) 

(4.23, x . (t „ )= ( 1 _^»5l^!l 

Notice that on the support of x T ■, we have |£ — 2r/| < ^ and hence, ygjj|£| < 2 1 77 1 < and 
|e-r?|<H + J|<^H. This yields 

3/2 /im\ 3 / 2 



r /200\ J/ " flOlY' \ ,3/2 

0->[(^J -'-(-W) M 3/2 >ct- 

In a similar way, we have on the support of x S that |£| 1/ ' 2 + |v7| 1//2 < |cL^__|. 

The rest of the argument is similar to what was done before. We rewrite the Duhamel term in 
(|4.4p as the sum of two terms: 

Y 2 m\° I \ X(^r ] )e i ^--^m(^r ] )f(s,r,)g(s^-r,)dr ] ds 



where x = X T or X = X S ■ F° r the term with we integrate by parts in time and recover terms 
similar to (|4.16b|) - (|4.16e[) . We estimate them exactly as in the case of the phase </>++. For the 
term with x ) we integrate by parts in rj, we get terms similar to (j4.17bp - (j4.17dj) . 
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d j r k Y 2 (d)ZA a f 

for some j, k less than respectively 7, K+3, which we fix from now on. First, since the commutators 
are controlled by (|4.11j) . it suffices to bound 

Y 2 {d)A a cPr k Y<f 



Since Y 2 (d), d, r commute with the group e ltA3/2 , we infer from (|4.1ip that 



Y' z (d)A a d j r k Sf 



\x- 



Since £ = S — §i$t, we can bound 



(4.24) 

Thus it suffices to prove that 



Y 2 {d)A a d j r k Y,f 2 < Y 2 (d)A a d j r k Sf 



+ t 



Y 2 {d)A Q d j r k d t f 



< e + ll«llv + t 



Y 2 (d)A a d j r k d t f 



Y 2 (d)d j r k d t f 



< \\ u \\x 



Going back to f|4.5j) . / solves 



d t f = e UA3/2 T m (u,u) + e UA3/2 R. 



Applying Y 2 (d)A a d^ T k gives 



(4.25) Y 2 {d)A a d' J r«d t f = e itA3/2 Y 2 {d)A a d j r k T m {u,u) + e ltA "'"Y' 2 {d)A a d 3 r K R. 

The remainder term R is easily dealt with, thus we discard it. 

Next, we can compute Y 2 (d)A a &> r k T m (u, u) by applying Leibniz' rule. This manipulation 
produces a great number of terms, but for simplicity we keep only two representative ones, namely 
(assuming for simplicity j, k even) 

Y 2 (d)A a T m (d j r k u,u) and Y 2 {d)A a T m {d j/2 r k / 2 u , d j/2 r k / 2 u). 

Applying Corollary 18.11 gives 



itA 3 / 2 



Y-'(d)A a T m (d j r k u,u) 



< 



A a d j r k % 



H 2 



A 3 /\ 



where the last inequality follows from (14. lip . Similarly, we find 



Y 2 (d)A a T rn {d j l 2 r k ' 2 u ) d j / 2 r k / 2 u) 



Coming back to (|4.25[> . this yields 



< 



A a dH 2 r k i\ 



H 2 



< \m\x 

w 2 <™ ~ t 



^3/4^/2 pk/2 t 



< \Mx 



Y' z {d)A a d j r k d t f 



< \\ u \\x 



thus giving the desired estimate. 

4.10. Proof of Q4.13P - It follows from ([4TT2]) , and Proposition \5l 
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4.11. Estimates for the remainder term R. Recall that we wrote in (|4.5p the equation under 
the form 

iA i/2 u + T m (u,u) + R 

with 



R 



d t u 

A 1 / 2 G{h)^ + i ( A h h + 



1 



(G(h)i/> + dh ■ Vi>y 



2(1 + 19/112) 

where [-]3 means the terms of order 3 and higher in the expansion of the expression between brackets. 
In proving the estimates (j4.8j) to (|4.13p , we always skipped the term R to focus on the more difficult 
bilinear term T m (u,u). We now come back to the term R and show how it can be estimated. We 
will prove the a priori estimate 



/' 

Jo 



Y(dyA a &r k R 2 ds < \\u\\ 6 x . 



for \j\ < 9 and \k\ < K + 8. We fix such a k from now on. Instead of treating all the terms which 
appear in the above expression of R, we retain one, [yl 1//2 G(/i)V'] „; this will alleviate the notations 
and the other terms can be treated similarly. Thus we want to show that 



Y(dfA a d j r k A 1 / 2 [G(/i)# 



ds< 



\x- 



Commuting A 1 / 2 and r k , interpolating, and applying Cauchy-Schwarz, we see that the above can 
be bounded up to lower order terms by terms of the form 



r k v [G{h)^\, 



ds 



vr k d j [G(/o# 



1/2 



ds 



1-6 



with < 9 < 1. 



Since the real problem is low frequencies, we focus on the first factor above. Expanding first 
G(h) into multilinear operators as follows from Proposition 17.11 and then applying Leibniz rule as 
described in Lemma 17.21 gives that 



r k cP [G(/#] 3 = r k & M n(h, ..-M) 



n>2 



is a sum of terms of the type 



^2 M n (d ji r h h,..., d jn r in h, d jn+1 r in+1 if) . 

ra>2 iiH |-in+l=fc 

with |u| + • • • + \i n +i\ < K + 8 and |ji| + • • • + |jn=i| < 9. Using the bound given in Proposition 17. 11 
each of these terms can be bounded by 



I \\M n (d ji r h h, d jn r in h, d jn+i r in + 1 ^) \\ ds 

Jo 



< c n 



< c:\\u\\ n x +1 



|^' i+1 r ii /i|| ...I 

Moo I 

Y(d)A 3 ^d ji r h u 
1 1 



d^+ v r ln h\\ a J "+ i+i r^ +i v L ds 

Moo M r M2 

Y(d)A 3 ^d jn r in 



V0 ,n+i r' w + 1 u|| 2 ds 



1 



TZsr(s) ds 



Summing over n gives the desired result since \\u\\x is small enough. 
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5. Dispersive estimate for the linear group 
Proposition 5.1. For any (3 6 [0, 



(5.1) 

In particular, 
(5.2) 



oo — ' — J 1 1 

j=0 fc=0 



1 3 



L 2 



1 3 



<7EE|M^~ 1/2Mfc / 



j=0 A;=0 



L- 



Remark 5.2. These estimates are probably far from being optimal, but, omitting Y(d), they are 
at the scaling of the equation. They should be understood as giving the right decay if f is smooth, 
and has a behaviour at oo in space (that is, at zero in Fourier) like an inverse power. 

A similar estimate was proved by Klainerman |23| . see also Hormander |22j for the linear group 
of the wave equation, and by Wu [M] for the group e ltAl/2 . These authors argued in physical space, 
whereas we rely here on a Fourier space approach. 

After some preliminary steps, the proof reduces the problem to an oscillatory integral, which is 
then estimated. 

5.1. Preliminary steps. 
5.1.1. Why (E2P follows from (El). Split 



e itA 3 / 2 j 



as follows, for jq < 0: 



/ 



< 



P<^f 



+ 



P>^ tA3/2 f 



Bound the first piece directly, using successively the Hausdorff- Young and the Cauchy-Schwarz 
inequalities, 



P<h e 



f 



< 







2-J° 



< 







230 



Id 



1/2-/3 



< 2 j °( 



A P-l/2f 



and the second one using ()5.2 

3 



P>_^f 



~ \ E \ P >M Y ( d ) A ~ 1,2 ^ k f\ 2 < \y{d)A^ 1 ' 2 T 1 Q h f 



k=0 



k=0 



Optimizing over jo gives the desired result. 



5.1.2. A Hardy-type estimate. We denote the polar coordinates in Fourier space by (p,0). A func- 
tion f(p) is understood to depend only on the radial variable. We will denote 



Mp) 



de_f J p 
P 



)V2-e if p < 1 
V2+6 if«>l. 



Then one can estimate by Cauchy Schwarz 

roo ^1/2 



\f(p)\ = 

We record the result 
(5.3) 



d p f(a) da 



< 



J) 



a p /((7) a 2 ^J < —wp^fw, 



\f{p)\<^\\ P 1/2 d P 7\\ L , {W) . 
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By a trivial modification of the above one can prove the slightly more precise bound 



(5.4) 



l/MI < ^ 



P(p)<VIIl2(r2)- 



5.1.3. Estimates on Bessel functions. Recall that the m-th Bessel function is given by 



3 i(s cos(8)+m0) 



It is well-known (see Stein [28], Chapter VIII) that 
J m (s)=e" s r(s)+e-"g?( a ) with 



< C(m,k) 



In order to bound the constant C(m, k), we will need the following lemma 

Lemma 5.3. Assume that f is a fixed smooth function such that /(0) = 0, /'(0) = 0, f'(x) ^ if 
x^O, and /"(0) = 1. Let 

I(s) := [ u(6)e isf ^ d6 



ifk = 0,1. 



where u is a smooth function with support in [—1,1]. Then 

k 



< 1 

~ ( s )l/2+fc 



l^lloo ~t~ ||oo ||^ 1 1 OO I 



In the definition of J m , the phase cos(#) has two stationary points, and tt, both of which are 
non degenerate. Successively restricting attention to one of them by an appropriate cut-off, we see 
that the above lemma implies 

Corollary 5.1. The constant C(m,0) and C(m, 1) can be bounded by 

C(m, 0) + C(m, 1) < m 2 . 

Proof. Choosing a compactly supported, smooth cut-off function, equal to one in a neighbourhood 
of zero, write 

I(s)= /(u(0) + 6u'(O)) X (0)e isfW d9 + [ [u{9) - (u(0) + 6u'(Q))x(0)] e isf(e) d9 := h{s) + I 2 {s). 



The term I\{s) is easily estimated by the standard stationary phase lemma. We are thus left with 
l2(s). In order to estimate it, it suffices to show that, under the assumptions of the lemma, and if 
furthermore u(0) = u'(0) = 0, then 

k 

l^lloo ~\~ ||^ ||oo ~\~ ||^ lloo) 



1 



< 



(this estimate is even stronger than needed). We prove this inequality for k = 0, the case k = 1 
being similar. The idea is of course to integrate by parts, to get 

i fu(s)f"(s) 



u(6)e isf( - e) d9 = 
which can easily be estimated by 

I'tol < ^ (Hoc 

This concludes the proof of the lemma. 

5.2. The oscillatory integral point of view. 



$&<P m d9 - 

f'{s) s 



f'(s) 2 

+ \\U lloo + \\U Woo) . 



□ 
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5.2.1. Reduction of the proposition to an oscillatory integral estimate. Expand the function / in 
the statement of the proposition in spherical (or rather circular) harmonics: 

m = E Tm{py mQ . 

meZ 

Suppose x, a point in physical space, has the polar coordinates (r,6o). Then 

J(rpco S (9)+me) e it P 3 / 2 f^^ pdpdd = J2e im9 ° f J m {rp)e itp3/2 J m (p)pdp 



e itA3/2 f 



( x ) = E e 

m 



where J m is the m-th Bessel function. By the results of Section 15.1.31 the problem reduces to 
estimating an integral of the type 



Atp 3/2 -rp) g{r p)f(p)pdp 



with 
(5.5) 



< 1 

~ ( S )l/2+fc 



for k = 0, 1 



(the other integral, with phase tp 3 ^ 2 + rp being easier, we skip it). Denoting 

de f r 3/2 



R 



the integral becomes 



and 4>{p) = p 3/z - Rp, 



Our aim will be to prove the 

Claim 5.1. If g satisfies \5.5\) . then with an implicit constant independent of R, 



^g(Rpt)f(p)pdp 



A{p)d p f{p) 



L 2 



5.2.2. How Proposition \5. 1\ follows from Claim [X71 The proposition follows easily from the claim. 
Indeed, using in addition the results of Section 15.1.31 it implies that 



< 

oo I 



(5.6) 



lj2 m2 \\ A (p)9 P up) 

m 

~ \ \J2 m6 \\ A (p) d pf™(p) 

\ m 



1/2 



L 2 



j=0 k=0 

5.2.3. Estimates on <j). Estimates on eft will be needed, we will record them below. First notice that 



b' vanishes for p = |i? 2 . This motivates the introduction of the new variable 



h = p — -R . 



We start with estimates on <//: 

4>'{p) ~ < 



ifp< 2 -£ 



f R 

h 

R Xi 9 

y/p if p > WR 2 . 



if lR 2 < P < 10i? 2 
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Estimates on d p ^ will also be necessary: 



lt P< — 



R 
Tr. 



if lR 2 < p < 10R 2 
if p > WR 2 . 



(5.7) 



Aside from direct estimates, the basic tool will be integration by parts using the identity: 

1 



Q e it<f>{p) _ it<j>(p) 
itcP'(p) P 



5.3. Proof of the claim (|5.1|) . the case p < -J^. We assume from now on for simplicity that 



A(p)d p f(p) 



1. 



We estimate here 



MWg(Rpt)f(p)pdp. 



5.3.1. The case R 3 > \. Using (15. 5p . (15.3p . and the assumption R 3 > j, it follows easily that 



i 

Ri 



e u ^g(Rpt)f(p)pdp 



i 

Ri 



5.3.2. The case R 3 < yLr. Split then 



loot • 

Rt e^(p)) g (Rpt)f(p)pdp= / t23 + * 



I- 



def 



I + 11 



Estimate for I The term / can be estimated with the help of (15, 5p . (15. 3ft : 



1 I\ < 



i 

t 2/3 



y/p~dp< 7' 



Estimate for 1/ Integration by parts with the help of (15.70 gives 



(5.8a) 
(5.8b) 

(5.8c) 
(5.8d) 
(5.8e) 
(5.8f) 



II 



it<f/ (t~ 2 / 3 ) 



g(Rtt- 2 / 3 )f(t-^) t 



2/3U-2/3 



J Rt 



r z lt ^ Rt 9(Rtp)Kp)pdp 

273 



" e u ^g(Rtp)d p f(p)pdp 

t 2/3 
1 



P73 

1 

Hi 
1 

2/3 



-g(Rtp)f(p)dp 



e it %( — )g(Rtp)f(p)pdp. 
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Observe that on the integration domain of 77, <p'(p) ~ y/p, d p y-^pjj^j ~ ^§721 an d g,g' ~ 1. Using 
in addition (|5.5p . f|5.3j) . (I5.4P and the Cauchy-Schwarz inequality, this leads to the estimates: 



— ~ t VF273 VF273 i 

/"St 1 1 /"Ht 1 

1^1 </ —Rt—pdp = R dp<- 
(5.9) ^ * 

1 1 / r \ 1/2 / /■ \ 1/2 -. 

ia < — P 9 p /(p)dp < - (j P A 2 d p f(p) 2 dpj \j a- 2 dp) < - 

(the term fl5.8e[) can be treated like (I5.8dp . thus we skipped it). 

5.4. Proof of the claim (15. ID . the case p > -j^. We estimate here 

poo 

e u ^g(Rpt)f(p)pdp. 

1 

m 

5.4.1. The case R 3 > \. Split the integral as follows 

00 _ f -R 2 /10 f .±B 2 ^Jf ,4^2 + ^/1 f . W B 2 r oc 

^(Mp)), ' 



e^)g(Rpt)f(p)pdp= + '+/ ' +/ _ + 



d = f I + II + III + IV + V. 



Estimate for J Integration by parts with the help of (|5.7p gives 
(5-lOb) + e M«_l f^W*) * 

/■- 1 

(5.10c) - / " e^—Rtg'{Rtp)f{p)pdp 

/■— 1 

(5.10d) - / 10 e ^ g(R t p)d p f(p)pdp 

r— 

(5.10e) - / 10 ^— g (Rtp)f(p)dp 

(5.10f) - e^dp (J^J g(Rtp)f(p)pdp. 
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Observe that on the integration domain of II, 4>'{p) ~ R, d p y-^/jj^J ~ W^/p ' 9^ ~ "7s' 5'( s ) ~ 
}m- Using in addition (|5.5p . (|5.3p . (|5.4p and the Cauchy-Schwarz inequality, this leads to the 



^372 
estimates: 



r— 

/"To 111 11 

Estimate for J/ Integration by parts with the help of (|5.7p gives 
(5.12a) 



J = -e 
(5.12b) 



M4) 1 



^' ( £ 



R 2 \ \ 10 / V 10 / m 



+ e 



it<j>[ 



(5.12c) 



r 



(5.12d) 



(5.12e) 



(5.12f) 



r 



ill 
10 



9 Rt 



ip2_. R 



it<f>> (|i? 2 - yffj 



e^^-Rtg'(Rtp)f(p)pdp 



-R 1 



-R 2 



4 r>2 R 



4 e u ^g(Rtp)d p f(p)pdp 

vtqy 



i 7?2_. R 



e^^-g(Rtp)f{p)dp 
it(f)' 



4 r>2 R 



4 i^ 7 )9(Rtp)f(p)pdp. 
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Observe that on the integration domain of II, <j)'(p) ~ |j, d p \jpjj^J ~ jp, g(s) ~ g'(s) ~ ^72- 
Using in addition (|5.5p . ()5.3|) . (j5.4j) and the Cauchy-Schwarz inequality, this leads to the estimates: 



,, .. ,1R 11 1 

m ~i7mvmn R =t 



(5-13) [r / \ 1/2 

,, .. ^ rvT i? 1 ,„ - Nl „ ^ 1 1 ( r°° d/i\ ' 1 1 1 



90 



i? 1 1 ^ VR f 00 dh ^ 1 



(notice that (|5. 12a|) was already estimated for /, and that (|5.12e[) can be estimated like (|5.12f|) . 
Estimate for III This term can be estimated directly, using simply (|5.5p and (j5.3|) : 



/•l^+Vf 1 1 1 



Estimate for IV Identical to the estimate for III 

Estimate for V Integration by parts with the help of (|5.7[) gives 

(5.14a) V = — it ^^ R2) 9 (RtWR 2 ) T(WR 2 ) 10i? 2 

(5.14b) - / e ^—Rtg'{Rtp)f{p)pd P 

f'°° 1 ~ 

(5.14c) - / e ^—g(Rtp)d p f(p)pdp 

J 10R 2 lt<pt 

(5.14d) - / e*» g(Rtp)f(p)dp 

J 10R 2 tt(// 

(5.14e) - J™ e^d p Q^j g(Rtp)f(p)pdp. 

Observe that on the integration domain of V, 4>'{p) ~ y/p, d p (^j) ~ ^72, g(s) ~ ^j, ff'(s) ~ ^72 . 
Using in addition (|5.5p . (|5.3p . (j5.4j) and the Cauchy-Schwarz inequality, this leads to the estimates: 
f 00 1 11 111 

iioi? 2 {tRpyi 1 y/p t s/tR? t 

r°° 1 1 ^ 1 / r°° dn\ 1 / 2 1 1 1 

(5.15) M< * |3 p /(p)|p d p<-i= / ^ 



,. ., „ f°° 1 1 1 1 ,^11 .1 

M < / 7 -^-==— P d P < 7- 



< - 



! WR 2 t p^VRtpVP ~ t v® ~ t 

(the term (|5.14ap was already estimated, and ()5.14dp can be estimated like (|5.14c|) ) . 
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5.4.2. The case R 3 < Integration by parts with the help of (15. 7p gives 
(5.16a, £<*™*MT(P)P*> = " ! *'^7TJS («£) T (±) ± 



nOO i 

it<j> 1 

US' 

Rt 

OO 1 



(5.16b) - / e lt f—Rtg'(Rtp)f(p)pdp 



(5.16c) - / i ^^-g(Rtp)d p f(p)pdp 

Tit 

oo y 




(5.16d) - / e u 't>^-g(Rtp)f(p)dp 

(5.16e) - j2 e^d p (J-^ g{Rtp)f{p)pdp 

Observe that on the domain of the above integral, 4>' ~ y/p and d p (^?jp)J ~ ^72 • Using in 
addition (|5.5|) . (|5.3|) . (|5.4p and the Cauchy-Schwarz inequality, this leads to the estimates: 

, , f°° 1 11 1 

f 00 111 1 

I WIS 4 Timvm7P pd "~l 

(the term (|5.16dp can be estimated similarly to (|5. 16cj) . thus we skipped it). 

6. Traces and elliptic estimates 

In this appendix we present bounds on harmonic extensions and traces of functions defined in a 
domain f Cl 3 . We assume that 3> is bounded by the graph of a smooth function h defined on 
R 2 such that 

(Al) \\A^+ a h\\ H2a _i_ a < 00, \\dh\\ Hs+ s « 1 for S »l, 0<a<i. 

This implies that h £ L p for some p < 00. For a function it defined on 

9 = {(x,z); x E R 2 ,z < h(x)}, 

we denote by u\,(x) = u(x,h(x)) the trace of u on the boundary S§ = dS 1 , and for a function 
</? defined on SB we denote by hs harmonic extension to Finally we denote by Aq g the 
solution to the Dirichlet problem on @ with boundary value 

Ait = g on S>, Uf, = 0. 

The surface S§ is a Riemannian submanifold of M 3 with covariant derivative given by 

(6.1) V e w = (V e w) T = V e w - (V e io • N)N = V e io + (10 • V e N)N, e,w £ TSS 
where N = |v(^-fe(a))| ' * s ^ ne um ^ norma l to In terms of coordinates {x l ,x 2 ) € M 2 

(6.2) ^ = {R 2 ; 9ij = 5ij + dihdjh, V i = d i + C i }. 
where Cj are the Christoffel symbols matrices. 
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Recall that Sobolev spaces on M? or R 3 are given by the norm 



\ u \\H k (R n ) 



{dy u 



L 2 (R") 

We will now define Sobolev spaces on T> and B. First, the function spaces H k {3>) are defined for 
any k as restrictions: 

llsllff fe (x>) =/ mf {llCll^fe^s) = #j . 

Notice that this definition is equivalent for k an integer to a more direct one, as can be seen using 
Lemma 16.21 below: 

Y H^sIIl^b) + \\dh\\ Hk - 1{R2) \\du\\ LOO(M?) ~ \\g\\ H k(T>) + II dhWgh-iQp) WduW^^ 

\s\<k 

For k <E N and < k < 2s, H k (m) is given by the following norm 

k 

ll ll \\T\i ll 

\\ u \\H k {SS) = / , \\ v u \\L^(m- 

By dSH, for 3 < k < 2s 

\\ u \\h*(R2) + ll^lliffe-l(R2) l|3w|| £0 o (K 2) ~ |M|ffk(B) + ll^ll^-i(B) ll^lli-(B) 

The next lemma shows that in the L 2 (33) norm, the Dirichlet to Neumann operator is the same as 
one derivative. 

Lemma 6.1. Given a harmonic function <j) on Si then the following are equivalent 

\\(9»<l>)b\\u> ~ \\(di<P)b\\ L 2 + \\(9 2 (f>)b\\L 2 ~ W&h* ~ I|V06||l2 ~ ||AMl 2 - 
Proof. Since is harmonic, then 

= J A(J)d z <j) = J N ■ V<J)d z <j) - ^N 3 \V(f>\ 2 

W 3 |0 Z 0| 2 = / N 3 [(di(l)) 2 + (5 2 </>) 2 ] - 2d^(iV 1 d 1 </> + N 2 d 2 4>). 



Since iV 1 and iV 2 are o(l) and iV 3 is 0(1) then we have the first equivalence. The second equivalence 
follows from the definition of M and from the relative sizes of N 1 , N 2 , and iV 3 

Na = n- (V0) 6 = iv 3 (a z 0) 6 + N^d^t, + iv 2 (d x2( />) fe . 

The third equivalence follows from 

di(j) b = (di4> + dihd z 4>) b . 

The last equivalence follows from using the Fourier transform and Plancherel's Theorem. □ 

The lemma above shows that for 1 < k < 2s we have 
(6-3) IMIjy*(0) ~ \W\\m-^(,gg) + ||^ fc-1 «||La(«)- 

In particular on H l ({%) we have the equivalent norms 

IMIl 2 (r 2 ) + II^IIl 2 (r 2 ) ~ IMU a (#) + II^IIl 2 (^) 

and by interpolation we can define 11^(9$) with equivalent norms 

IMU 2 (R 2 ) + |I^ 5 V||L2(R2) ~ ||<^|| L 2(^) + \\N 2 if\\ L 2 . 
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The spaces H k+l / 2 {3§), for 1 < k < 2s - ±, can be defined by 

(6-4) \\(p\\ H k+i/2 {m) d = \\<p\\ H *(a) + \W^V k ip\\ L 2 m . 

To show that harmonic extensions gain 1/2 a derivative and that Aq 1 gains 2 derivatives, in an 
appropriate range of spaces, we introduce an H 2s+l / 2 coordinate system that maps M?_ — > & and 
the z = plane to S3. 

Lemma 6.2. 1. There exists an extension h — >■ h defined on M 3 such that 

||V/l||#2 s -l/2( R 3) < ||d/t||H2 s -l( R 2) llV/tH^+l^^S) < ||9/l||H«(R2) 

2. The map y = (x,z) y = (x, z) defined by 

x = x 

z = z + h(x, z). 

the z = plane to S3, and maps the lower half space Ml — > *2) . ^ is an H 2s+1 ^ 2 diffeomorphism 
on M 3 , since 

(6.5) H-D* - I|| H 2s-l/2 < OO, ||D* - I|| H s+7/2 < 1. 

Proof. Define h via its Fourier transform h (£ denoting the dual variable to (x%, X2) and n the dual 
variable to z) 

ktr]] ,= :c WT7^ h{0] c " 1 = J " {i + s 2 r+^ ds - 

Then h(x, 0) = h(x) and for 1 < a + b < 2s + 1/2 

(M) b 

\m b HM = icr 6 - 1 - — ^ l \ a+1/2 ko => \m a \v\ b h\\ L z < m a+b - 1/2 h\\^ 

M n2 



1 _i_ 1 m 

i + Ui 

By assumption (lAlj) and the definition of h 2. follows. □ 

The next proposition states that for < a < 2s, 1. g — > g o \s one to one and onto from 
to H a+ ^(R 3 _), 2. harmonic extension of / E H a (38) gains 1/2 derivative, and 3. VA^V 
is bounded on H a ~^(S#). 

Proposition 6.3. Under assumption (|A1|) and for a G N/2 and 3 < a < 2s we have: 

1. For g G H a+ ^(&) 

\\g\\ Ha+im < h O ^\\ H a + ^ m + \\dh\\ Ha - 1{m \\V(g O *)|| w[ a ]+1)0o(H 3^ ) 

2. For f G fP(^) 

llV/^ll^i ^ < \\f\\ H » (3g) + ll^ll^-i^HV/ll^]^; 

3. For g G # a_1 / 2 (f^) 

HVAo^pll^i <IMI^-i ^ + I|S/»||h-i W ||VA^ 1 V & || w14]+1 , 00( ^ ) 
^Nl^-^ + Hd /l lltf-i(^)ll V 2ll H [f]+ 2 (#) 
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Proof of 1. For a + 1/2 an integer the inequalities follow immediately from the chain rule. For a 
an integer, extend g to a function g e £ fC a+1 / 2 (R 3 ) with ||g e ||_fl-a+i/2( K 3) < 2||g||^a+i/2(^) and use 
the equivalent norm on fC a+ 2(R 3 ) 

f f \D a ^{x) - D a <t){y)\ 2 , , 
H 2 J j \x — y\ 

and (16.51) to conclude 1. 



Proof of 2. and 3. For g £ H a ^^{&), solve AC/ = Vg on ^ with Dirichlet data /. Change variables 
using ^ to flatten the domain, and let <fi = U o ^, an d j = j o $; then DC/ = ((/ + A)D(f>) o if -1 
where the matrix A depends on <9/i, and thus satisfies 

(6.6) Il^ll.ff2s-1/2( R 3) < OO, 11^11^+5/2^3 ) < 1. 

Consequently (j) satisfies 

trace[(L> + AD) 2 ](j) = A<J) + 2 div(5Z)^) + c- D<j) = Vg + AVg i£l 3 

0(x,O) =^b(x) =/o*(x), 

where c depends on d 2 h, and i? is a linear combination of A and A 2 . By letting A" 1 be the inverse 
of A on Rl , we write 

D<f> = D(<p b ) H + DA^Vg + DA'^BVg - c D<j>) - 2VA; 1 div BD<fi. 

By standard elliptic estimates on R 3 , and Sobolev embedding we have 

WWWvfpL) £ IMfrhrn^ + Mvpl) + \\BD0\\ L2m + \\BVg\\ « + ||c- £>0||. fi 



IIV^II^-l^ } < WMh'V?) + 11^11^-1^3) + Hffll^-l^J 

+ ll^ll^-a^) + \\BVg\\ L i m + He • ^ll^-| (R 3_) + He • m L ^ R3 _y 



By setting g = 0, or / = 0, and using (|6,6p . we conclude 2. or 5., for the stated ranges of a. □ 
The next proposition shows that 

g 

i=0 

Recall that di = di + hid z , and that for a vector field X = X ■ V on T^, the commutator [<-f ,A/] 
acting on cp : SB — > R is given by: 

(6.7) = N • VA 1 (2div((VX w )V^)) + (XN) ■ (X7 m ) - N • {DX n D m ) 
Thus by proposition 16.31 

(6.8) \\[d,Ar\<p\\i? ~ \\(D 2 h)D<p\\ L i 
Proposition 6.4. For any ip : SB — > R and 3 < q < 2s with q G No/2 

I|aavIIl 2 < ELo II^IU^) + ll^n^-x nv^Hicc 

II^Mtf < ELo II^IU^) + \\dh\\ m -i \\V<p\\ L ~ 
Proof. From [29] we have 

(6.10) - N 2 )y = kN{<p) + 2N • VA 1 dw{DN n Dip n ) - Af(N) • (tf(<p)N + V T ip). 



(6.9) 
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Thus the operator A^+A^ 2 is a first order operator with coefficients depending on second derivatives 
of h. Consequently if q is even (16. 9ft are immediate. If q is odd then the equation (16. 10h and Lemma 
16.11 imply (|6.9p . For q = k + \ , we use the commutator estimate f|6.8j) to write 

\\Ni<pf L2 = J J\f h+1 (fJ\f k tf = J D k ipND k <p + lower order terms. 

^From equation f|6.4[) we obtain the desired result. □ 

7. Properties of the Dirichlet-Neumann operator 

Recall that we denote M{h) or simply J\f for the Dirichlet-Neumann operator associated to the 
domain T> = {z < h(x)} with boundary SB = {z = h(x)}. 

7.1. L°° estimate. Next we bound M in L°°{SB) = L°°(IR 2 ) by allowing the loss of small powers 
of A = |V| for small and large frequencies. 

Proposition 7.1. There exists eo such that: if 

Halloo + l|V 3 /i||oo + \\dh\\ 2 + ||V 3 /i|| 2 < e , 
then for any (p defined on SB, and for < a < 1/2, 
(7.1) HA/VH^wi^ < ll^llws.oo + ||^>||l 



oo . 



Proof. Step 1: the double layer potential Let a = (x, z) £ f and 6 = (y, /i(y)) G Then can 
be represented by the double layer potential 

m (a) = j n(b)N ■ VG(a - b)dS(b) = -l/2fi(b ) + J (p(b) - (J,(b ))N ■ VG(a - b)dS(b) 
m a 
where G(a — b) = j-|a — b]" 1 is the Newtonian potential and bo is an arbitrary point on SB. Define 
K(x, y) = y/l + \dh{y)\ 2 N(b) • VG(6 - b) for b = {x, h{x)) G SB and b = (y, h{y)) G In other 
words, 

= -V%) • (a - y) + - %) 
4vr(|x - y| 2 + (h(x) - h(y)) 2 )^ 
Then as a — > bo, it follows by standard singular integral calculations |14j that 

(7.2) - l/2/i(x) + f n(y)K(x, y)dy = <p(x). 



Step 2: estimating \i Notice that K(x,y) has the following properties: 

• / K(x,y)dy = 

• \K(x, y)\<e min , -^p 

• \K(x, y) - K(x', y)\ < eoj^rfpr if \v ~ x \ » \ x ~ x '\ 

• I \ K (x,y)\ dy < e . 

The first point above can be seen by a standard integration by parts; the second and third points 
follow from H^/illoo < eo; for the fourth one use in addition that ||<9/i||2 < eo- 

Next we also denote K for the operator with kernel K(x,y). By points one to three above, it is 
standard to see that 

H-Kllca-wi" ~ e ° if < a < 1; 

and point four easily implies 
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Choosing eo small enough, it is possible to solve (|7.2p by Neumann series to obtain 

IIHk 00 ~ IMk°° an d Wfj-Wcc < \\<p\\c<* 

To obtain estimates on derivatives of //, write 

-l/2fi(x) + / (fj,(y) - n(x))K(x, y) dy = ip(x). 



and differentiate with respect to x to get 

-l/2d x fi(x) + / (/i(y) - fi(x))d x K(x,y)dy = d x tp(x). 



Denoting J(x, y) = (9 X + d y )K(x, y), an integration by parts gives (writing J for the operator with 
kernel J(x, y)) 

(7.3) - l/2d xf i(x) - Kd xf i + Jfi = d x ip. 

Using that ||V 2 /i||oo + ||V 2 /i|| 2 < eo, one checks that J(x,y) enjoys the properties of K(x,y) listed 
above. Thus if < a < 1, 

(7.4) ||</||l°°->l o ° + || J\\(i a -+ca < e o- 
It is furthermore easy to see that 

ll<^llc 1 -' T ->-L 00 + II^IIci-ct_j.loo+ ^$ CO- 

Now equation (17. 3j) implies 

d x fi(x) = y^ + KJ (Jfi-d x <p), 
from which follows, in conjunction with the various bound on K and J, if < a < 1, 

Halloo ~ Halloo + Halloo < Ibllci-- + Halloo 
II^bMHci-o- < HJ/illci-o- + H^cVlIc 1 — " ~ ll^lld 1 -- + II^Vlloo- 
Taking one more derivative and arguing similarly gives 

II^mIIoo < WdlvWoc + IMIci— • 

Step 3: estimating N<p Fix a point 6o € S3 and use normal coordinates in a neighborhood of SB to 
restrict a near the boundary to the line a = bo + vN(bo). Thus 

jV(6 ) • Vp«(a) = J (jt(b) - fi(b ))D 2 G(N(b),N(b ))(a - b)dS(b). 

ss 

For 1 1> — £>o I large and v small \D 2 G(a — b)\ < \b — 6o|~ 3 an d thus the above integral can be bounded 
by ||//|| c i 1 _ CT . For \b — bo\ small we write 

N{b Q ) = 0(6, b )N(b) + 7(6, b )r where r G T b S. 

The term involving r is integrable due to the vanishing of 7(60,60)- By repeating the argument 
that led to inequality (|7.4|) we obtain 



(H(b) - ^bo))D'G(N(b) iT r)(bo -b)dS(b)\ < M d i-« < \\<p\\fr- 

Split the remaining term 

1=1 0u(6) - v(b ))D 2 G(N(b), N{b))(uN(b ) + 6 - b)dS(b) 
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with the help of a cutoff function % which localizes smoothly away from the ball of center and 
radius one: 

I = J (x(b - b ) + [1 - x(b - b )]) . . . dS(b) = 11 + III. 
The term III is easily dealt with; as for II, it is hypersingular as v — > 0. Use the identity 

= AG = A S G + kN • VG + D 2 G(N, N) 
for v < to re-express as 

11= [v ( x (b - bo)(ji(b) - fi(bo))) VG(uN(b ) + b -b) 



- X (b - b )(Kb) ~ fi(b ))KN(b) ■ VG(vN(b ) + b -b) dS(b). 
Letting now v — > 0, this can be bounded as before 

ll-HI ^ ll^llci + ll/^llci— ^ ll^llci + IMIdn-*- 

Thus 

\\N<p\\l°°(S) < WfyWc 1 + IMIc 1 — ■ 
By repeating the above argument after applying tangential derivatives to N(bo) ■ V</> we obtain 

\\M^\\ W 2,oo {s) < \\d<p\\ca + ll^llc-i— 

This proves inequality (|7.ip . □ 

7.2. Expansion into multilinear L 2 bounded operators. 

Proposition 7.1. The operator M(K) can be expanded into a series of multilinear operators 
N{h) = Y J M n (h,...,h,-) i.e. Af(h)f = Y,M n (h,...,h,f) 

n>0 n>0 

where the operators M n are symmetric, n — linear in h, and satisfy the estimate 
(7-5) \\M n (h,...,h,f)\\ 2 <C?\\dh\\^ ||V/|| 2 . 

Proof. Though we could not find this explicit statement, this result is classical in harmonic analysis. 
The idea is to expand the Dirichlet-Neumann operator into elementary operators akin to Calderon's 
commutators. 

Step 1: The single-layer potential Let ip be the solution of the Dirichlet problem with boundary 
data /: 

= in T> and ip(x,h(x)) = f(x) on SB. 
It can be represented via a single layer potential by 

i>(x,z)= [ ^ -^dy 

Jr* (\x - y\ 2 + \h(y) - z\ 2 ) 1/2 



(this is not quite the standard definition, where the weig ht ^l + \dh\ 2 would appear, reflecting the 
volume element for SB). In particular, 

p(y) 

(\x-y\* + \h(x)-h(y)\*y 

It is well-known (see Folland [13]) that 

M{h)tp{x) = lim^,,^) N x ■ Vip(z, h(x)) 
(7.7) 1 p(x) f h(x) - h(y) + 8h(x) -(x-y) 



2 Vl + |5/i| 2 (x) V (\ x -y\ 2 + \h(x)- h(y)\ 2 f /2 ' 
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Step 2: expanding p Start by expanding (|7.6p in h: 

where |/3 n | < 1. Apply A to the above, noting that A = -j • D = X^'=i 2 ^i^i- This gives 



n=l j=l,2 



X — y| 2n + 1 



2„a J ./ tW -(2» + i)'f-^" 1 W-y 

|x — y| | a; — y| 



def 



p^ + Y^lPnP] (X). 
n=l 

Inverting the above by Neumann series gives 



(7.8) p(x) = £(-!)< 



fc=0 



n=0 



/!/(*). 



Step 3: expanding M(h)f The right-hand side of (j7.7j) can be expanded in h to give 
(7.9) 

1 \- h(x) - h{y) + g%) -(x-y) \h{x) - h{y)\ 2n 

[N {h)ip\ (x) = -—======== + > a n : 13 : p(y dy, 

2 yq + |Mp(x) ^ |x-y| 2n 

tfe/ 1 p{x) 



oo 



2^l + \dh\ 2 (x) 
where \a n \ < n. 

Step 4: boundedness of the elementary operators and conclusion The elementary operators P n and 
Q n appearing respectively in (|7.8p and (|7.9p can be estimated by the following result of Coifman, 
Macintosh, and Meyer [9] 0: the operator on M 2 with kernel 

- fel(y)) • • • (M^) ~ My)) 

| a; — y| n + 2 

(and n odd) has a norm on L 2 which is bounded by 

ll || oo • • • || Ofl n || oo • 

Using this estimate in conjunction with (j7.8|) and (|7.9p gives the desired conclusion. □ 
7.3. Symmetries. 

Lemma 7.2. The Dirichlet- Neumann operator is invariant by the following symmetries: 

• Translation: G{h(- + 5)) [/(• + 5)} = [G(h)f] (• + 5) for 5 e E 2 . 

• Rotation: G(h(Rg-)) [f(Rg-)] = [G(h)f] (Re-) with 9 € K, where Rq is the rotation by an 
angle 9. 

• Dilation: G (jh(X-)) [/(A-)] = A [G(h)f] (A-) toitt A > 



c(i + n 4 )||a/i 1 || 00 ...||a/ lr 



^This result corresponds to Theorem III in [9], transferee! on R 2 by the method of rotations. 
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Consider the expansion of Proposition 17.11 Applying for instance the rotation symmetry, and 
using in conjunction the above lemma gives 

[G(h)f] (R e -) = £ M n (h(R e -), h(R -), f(R e -)) 

n>0 

Differentiating in the continuous parameter [9 in the above example, or A or 5 if the other symme- 
tries are used), we obtain that 

r [G(h)f] = nM n (rh, ...,hj) + M n (h, ...,h,rf) 

n>0 

for r = V, i?; and 

(7.10) r [G(h)f] = Y, nM n (rh, ...,h,f)- (n + l)M n (h, ...,h,Tf) + G(h)f. 

n>0 

for r = x- V. This last formula remains valid for S = |i(9j+x-V if the functions under consideration 
depend on time. 

Analogous formulas can of course be obtained if more than one vector field P is applied. 

8. Tools from linear and multilinear harmonic analysis 

8.1. Littlewood-Paley theory. Consider a function supported in the annulus C(0, |, |) such 
that 



j 6 

Also define 



for^O, £*(|)= L 



j<0 

Define then the Fourier multipliers 

«Hv) —(f) — <§ 

and similarly P<j, P>j- This gives a homogeneous and an inhomogeneous decomposition of the 
identity (for instance, in L 2 ) 

Y p j = M and p <o + Yl p j = Id ■ 

All these operators are bounded on L p spaces: 

ifl<P<oo, ||P,-/II P < ll/ll P , ||P<y/|| P <||/||p and ||P>,-/|| P < ||/|| p . 
Furthermore, for Pjf, taking a derivative is essentially equivalent to multiplying by 2- 7 : 

if l<p<ooandaGM, \\A a Pjf\\ p ~ 2^'||P i /|| p 
if 1 < P < oo and £ G Z, \\V e Pjf\\ P ~ 2^'||P i /|| p . 
Also, we recall Bernstein's lemma: ifl<<?<p<oo, 

(8.2) < 2 2j (!-i) \\ Pj f\\ and ||P<,/|| p < 2 2j (^) \\ P<j f\\ . 
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8.2. Boundedness of bilinear operators. Recall the definition of the bilinear operator with 
symbol m: 

T m (f, g)(x) = f [ e**f(ri)g(Z - 7 ? )m(£, 77) d^drj. 



These operators are called pseudo-products and were introduced by Coifman and Meyer [10] . 

Definition 8.1. Let Jvl^^ 1 ^ 2 ^ 3 denote the set of bilinear symbols m(£, 77) such that 

• m is homogeneous of degree (3: m(\£,\rj) = A^m(£, 77). 

• m is smooth away from {£ = 0} U {rj = 0} U {£ — rj = 0}. 

• If |£| ^ \v\ ~ 1; m (£> V) can be written under the form \£\ Cl A T]\ , where A is smooth 
in its arguments. 

• If\v\ ^ |£l ~ 1j m {i) r l) can be written under the form \i]\ C2 A (\rj\, j^,(^J , where A is smooth 
in its arguments. 

• If |£ ~ v\ ^ ~ I; ^(C) 7 ?) can be written under the form |£ — ij\ cs A M£ — r]\, , r] \ , 
where A is smooth in its arguments. 

Proposition 8.2. Letij)(£,r)) be a smooth function, supported on an annulus, and letm E .A4^> C1 ' C2 > C3 , 
with ci,C2,C3 > 0. Finally set 

Then 

(8-3) \\T»{f,g)\\ p <tf j \\f\\ q \\g\\r 

if l <P,q,r < oo and | + ± = ±. 

Proof. By scaling, it suffices to treat the case j = 0, (3 = 0. Consider thus the symbol tpfx. It is 
supported on a compact set. First consider the region inside this compact set where none of £, 77, 
£ — 7/ vanish. Then ^/x is smooth, and the proposition is clear. We are left with the three regions 
where \tj\ <C |£|, |£| <C |r/|, and |£ — 77 1 <C |t7|. By a duality argument and symmetry, it suffices to 
treat one of these three cases, say \rj\ <C |£|. Then m can be written 

\v\ e *(\v\>fryt 

Expand A in Taylor series in its first argument, and in Fourier series in its second argument A £ S 1 
(we omit the necessary cut-off function to alleviate the notations) 

/ \ M 

Aw> e»0 = EE H^ e *^(f) + 

^ ' keze=o 

(notice that the smoothness of A entails fast decay of the Am in k, £). Taking M sufficiently large, 
Rm is sufficiently smooth for the theory of Coifman-Meyer to apply, and we are left with Yle=o- 
We might as well consider only the first summand, I = 0, the other being treated similarly. Thus 
we are left with 

EM^^)- 



Multiply this symbol by ip(£,,i]), and call the result /?(£, 77) = ip{£,,rj) zJfcez \v\ Ce M -^ko{C)- Then 
T P {f\g) = E T„{P m f,g) = E E^tt.^oCO^*** Pmf,Pig). 

m<0 m<0fceZ 
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Now by standard linear theory there exists C such that 



A c e ik ^P m f 



< 2 mc k c 



valid for any p in [l,oo]. Summing over m is possible since c > 0; and the fast decay of the Ako 
makes the sum over k converge. This finishes the proof. □ 

A simple consequence of the previous Proposition is the following 

Corollary 8.1. Let m G M^' Cl ' C2 ' C3 and assume that °"3; Q, r , Q, R satisfy 

11111 

ci > 0, a 2 < c 2 , cr 3 < c 3 , and - + - = — + — = -. 

q r Q H p 



Then for any k > 0, 



\T m (f,g)\\ 



LP ~ 



[A K + /T K ] A^-^g ^ + [A K + A~ K ] A^-^f 



LQ 



\9\\w"s,,R 



Proof. As in the proof of Proposition 18.21 it suffices to treat the case where m is supported in the 
region where \rj\ « |£|. Then n^fep £ ^^-^.ci^-o^.ca Keeping the notation tp defined above, 
and applying Proposition 18.21 one gets 



<2^-^'|||r ? rVIUI^I 



LP 



Summing over j gives the desired result. 



□ 
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